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EXECUTIVE  SUMMARY 


A  principal  component  for  the  desiqn  of  future  surveillance  systems 
is  a  detailed  description  of  rocket  plume  characteristics.  This  description 
of  plumes  is  primarily  built  upon  state-to-state  excitation  cross 
sections  for  certain  very  specific  collision  processes.  These  processes, 
predominantly  involving  atomic  oxygen  and  plume  exhaust  gases,  are  now 
believed  to  be  the  crucial  driving  mechanisms  for  the  enhanced  intensity 
of  infra-red  radiation  from  plumes  at  altitudes  above  100-200  kilometers. 

Such  c-oss  sections  are  requisite  input  to  sophisticated  hydrooynamic 
plume  modeling  codes,  which  require  kinetic  information  for  as  many  as 
one  hundred  different  microscopic  molecular  processes. 

For  many  of  these  kinetic  processes,  the  required  information 
is  not  available  from  experiment.  In  other  cases,  the  experimental  methods 
have  proven  expensive,  and  sometimes  unreliable.  !t  is  felt  that  the 
sophisticated  methods  of  modern  theoretical  chemistry  can  often  provide  a 
cost-effective  means  for  predicting  state-to-state  cross  sections  caused 
by  collisions  in  plumes.  Theoretical  calculations  also  assist  in  the 
interpretation  of  the  experimental  results. 

The  intent  of  this  project  is  to  use  current,  state-of-the-art 
ab  initio  theoretica1  chemistry  computations  to  predict  a  series  of  cross 
sections  for  vibrational  excitations  caused  by  0(JP)  atoms  colliding  with 
the  two  plume  soecies,  H^O  and  CO^.  These  studies  are  the  first  results 
fo1-  systems  of  this  complexity,  determined  at  the  degree  of  sophistication 
employed  in  this  investigation.  All  calculations  are  "first-principle". 


7 


meaning  they  are  built  solely  upon  the  laws  of  physics  and  quantum 
mechanics.  Only  a  knowledge  of  the  atoms  and  molecules  involved  in  the 
collision  is  required  information.  There  is  no  scaling  to  experimental 
values,  no  adjustable  oararneters ,  and  no  attempt  to  use  crude,  often 
structureless, model s  for  the  collisions.  Furthermore,  most  of  the 
results  reported  here  have  been  obtained  prior  to  any  experimental 
information. 

3  1 

In  the  collision  of  0(  P)  with  HgOt  ) ,  results  are  reported 
for  the  vibrational  excitation  cross  sections,  (001),  (100),  (010),  (011), 
(110),  (200),  (020),  and  (002),  at  a  collisional  velocity  range  above 

4  km/sec.  For  C02,  the  cross  sections  (001),  (100),  (010),  (020), (110), 
(030)  over  a  velocity  range  of  6-10  km/sec  are  given. 

The  report  describes  the  three  steps  in  the  calculation,  the 
ab  initio  computation  of  points  on  the  potential  energy  surface  for  the 
collision,  the  method  of  making  an  analytic  fit  to  the  computed  points, 
and  the  details  of  the  cross  section  determination  by  quasiclassical 
trajectory  methods.  Emphasis  is  placed  on  the  possible  errors  that  can 
occur,  and  the  justification  for  the  reliability  of  the  results.  Several 
recommendations  for  future  work  are  offered. 
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I .  INTRODUCTION 


Above  100-200  kilometers,  rocket  exhaust  plumes  show  a 

significantly  increased  level  of  in  , ared  emission  compared  to  what 

would  be  expected  based  upon  assumptions  of  local  thermodynamic 
n  2) 

equilibrium'  ’  '  A  large  and  diverse  data  base  obtained  by  field 

observations  and  laboratory  experiments,  assisted  by  gas-dynamical  models 

of  the  upper  atmosphere,  have  explained  this  increased  infrared 

(1  2) 

emission  by  a  series  of  basic  physical  mechanisms.  ’  '  One  of  these 
important  mechanisms  involves  collisions  among  atomic  oxygen,  the 
dominant  lower  ionosphere  species,  and  plume  exhaust  gases  such  as  CC^ 

H2O,  HF,  and  F  .1.  These  collisions  cause  an  increase  in  the  population 
of  excited,  infrared  active  vibrational-rotational  states  for  the 
molecules  ’n  the  exhaust  gas,  which  then  emit  radiation  in  the  infrared  . 

However,  in  spite  of  our  understanding  of  the  basic  physical 
mechanisms  underlying  the  observed  phenomena,  there  is  little 
detailed  signature  data  for  low-thrust/low-intensity  targets.  This  lack 
of  information  has  made  the  design  of  future  surveillance  systems  quite 
difficult/1 

As  an  alternative  to  a  knowledge  of  detailed  signature  data, 

high-altitude  target  prediction  computer  models  have  been  used  to 

determine  the  infrared  intensity  levels.  These  computational  models 

employ  fluid  mechanics  for  nonlocal  thermodynamic  equilibrium  environments, 

and  among  other  things,  must  incorporate  data  for  nonenui 1 1 brium 

(2  3) 

molecular  radiation  phenomena.  •’  Consequently,  these  hi qh-al ti tude 
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signature  codes  are  based  upon  detailed  information  about  the  inelastic 
energy  transfer  processes  that  occur  in  hyperthennal  collisions. 
Specifically,  it  is  necessary  to  know  the  probability  that  a  molecule 
will  undergo  a  particular  vibrational-rotational  transition  when  colliding 
with  atomic  oxygen.  These  probabilities  are  referred  to  as  "cross- 
sections,"  and  these  cross-sections  are  required  input  to  the  computer 
models  of  plume  radiation. 

In  the  absence  of  experimental  determinations  of  such  cross 
sections,  it  has  been  necessary  to  rely  upon  estimated  and  probably 
unrealistic  values  for  these  quantities.  This  is  the  most  serious 
limitation  of  the  high-altitude  plume  signature  models.  Consequently, 
a  substantial  effort  toward  resolving  this  problem  by  obtaining 
reliable  values  for  these  inelastic  cross  sections  is  of  high  priority. 

Cross  sections  may  be  obtained  in  principle  by  shock-tube  and 

(41 

crossed-molecular  beam  experiments.  '  However,  such  experiments  are 
difficult  to  perform  and  to  interpret.  It  is  also  possible  to  predict 
cross  sections  from  first-principle  theoretical  calculations. 

This  is  the  intent  of  the  present  project. 

The  contribution  of  accurate  theoretical  calculations 
to  the  problem  of  vibrational  excitation  cross  sections  is  three-fold: 

(1)  Strong  theoretical  support  is  required  to  assist  in 
interpreting  the  experimental  data 

(2)  Theoretical  calculations  will  often  be  the  only  way 
to  predict  cross  sections  not  presently  accessible 
to  experimental  techniques 
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(3)  The  reliability  of  molecular  calculations  has  pro¬ 
gressed  sufficiently  in  recent  years  that  one  can 
anticipate  that  a^  initio  theoretical  calcula¬ 
tions  can  predict  accurate  values  for  inelastic 
vibrational  cross  sections  in  many  cases. 

If  such  theoretical  calculations  are  found  to  be  accurate  to  within 
a  required  tolerance  (nominally  a  factor  of  2),  then  the  calculations 
may  become  the  most  cost  effective  approach  to  providing  the 
requisite  cross  sections. 

In  the  followinq,  we  report  on  the  accurate,  first-principle 
predictions  of  the  vibrational  excitation  cross  sections  for  the  reactions 

0(3P)  +  H90(ooo)  -  v ^  ->i v^)  +  0(3P)  (1) 


and 


0(3P)  +  C02(ooo)  -  C02(v1v2v3)  +  0(3P)  .  (2) 

These  calculations  employ  techniques  that  are  state-of-the  art  for 

inelastic  collisions  of  this  complexity.  Furthermore,  with  the  exception 

of  a  single  shock-tube  ,  _.>ult  for  a  cross  section  in  H20  near  threshold 

( 5 ) 

that  dates  back  to  1975,  all  results  reported  here  were  obtained  prior 
to  any  experimental  information  on  these  two  collision  systems.  As  such, 
this  effort  represents  something  of  a  milestone  in  the  application  of 
first-princ  pie  theoretical  calculations  to  the  problem  of 
immediate  importance  to  the  national  defense. 
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The  theoretical  approach  to  the  prediction  of  cross  sections 


nas  three  parts: 

(1)  Prediction  of  the  Potential  Energy  Surface  (PES) 
for  tne  collision. 

( T J  Analytic  fit  to  the  PES. 

(3)  Uuasiclassical  trajectory  calculations  to 
obtain  the  cross  sections. 

In  order  to  obtain  high  accuracy  in  the  predicted  cross  sections,  it  is 
necessary  to  maintain  a  high  standard  of  reliability  in  each  phase  of 
this  study.  Section  II  provides  a  synopsis  of  many-body  perturbation 
theory  (MbPT)|6,7^  which  is  the  ab  initio  method  used  to  generate  points 
on  the  PES.  This  is  followed  by  Section  III,  a  detailed  discussion  of 
the  0(JP)  +  H^0( 1  A-j )  cross  sections,  subdivided  into  four  components, 
tne  PES,  the  analytic  fit,  a  discussion  of  the  quasiclassical  trajectory 
method,  and  the  actual  determination  of  the  vibrational  excitation  cross- 

*3  1  + 

sections.  Section  IV  reports  on  the  0(  P)  +  r^)  cross  sections, 

also  divided  into  four  parts.  The  final  section,  V,  summarizes  the  project, 
evaluates  the  accuracy  of  the  calculations,  and  provides  suggestions 


for  further  work. 


II.  A  SYNOPSIS  OF  MANY-EODY  PERTURBATION  THEORY 


The  calculation  of  several  points  on  a  PES  requires  accurate 
and  efficient  ab  initio  methods  that  properly  include  the  effects  of 
electron  correlation.  This  work  employs  MBPT  built  upon  an  unrestricted 
Hartree-Fock  (UHF)  reference  function  for  this  purpose. ^  Although 
comparatively  new  in  applications  to  potential  energy  surfaces,  MBPT 
offers  a  number  of  features  that  recommend  its  development  for  this  type 
of  problem. 

These  features  include  "size-extensivity"  or  the  correct  scaling 
of  the  calculation  with  molecular  size.^’^  This  is  important  in  exten¬ 
sions  to  larger  molecules  and  in  predicting  the  correct  dissociation 
energies  when  comparing  a  molecule  and  its  fragments . A  truncated 
configuration  interaction  (Cl)  calculation  is  not  size-extensive,  although 
"full  Cl"  (including  all  excitations)  and  certain  other  models  like 
generalized  valence  bond  (GVB)  are.  Also,  unlike  Cl,  MBPT  permits  the 
inclusion  of  quadruple  and  higher  Cl  excitations  in  convenient,  tractable 
ways  which  also  serves  to  enhance  the  comparative  efficiency  of  MBPT 
over  Cl. 

There  are  also  disadvantages,  largely  due  to  the  present 
practical  limitation  to  single  reference  functions.  (The  multi  reference 
function  MBPT  theories  exist, tiut  have  not  yet  been  implemented 
in  general  purpose  codes.)  A  restriction  to  a  single  reference  function 
usually  requires  that  an  unrestricted  Hartree-Fock  (UHF)  function  be 
used  to  study  open-shell  systems.  This  has  the  effect  that  for  some 
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cases,  large  degrees  of  spin-contamination  can  remain  in  the  calculations 
[Ng( ^ -g)  a  case  1n  point.  where  UHF  is  used  to  obtain  correct  separation 
to  N(4S)  atoms^],  although  for  many  genuine  open-shells  this  is  not  a 
problem  [e.g.,  0(^P),  N(4S),  HC0(^A')][^  In  the  present  calculations, 
the  multiplicity  of  the  0{°P)  +  HgO^A-j)  and  0(^P)  +  CC^^ig)  surfaces 
a^e  monitored  and  found  to  always  be  better  than  3.01,  so  the  UHF  starting 
point  does  not  suffer  in  this  application  from  excessive  spin  contamination. 
Also,  only  the  lowest  state  of  a  given  symmetry  can  normally  be  obtained 
from  the  single- reference  function  approach.  We  will  return  to  this  last 
point  in  the  discussion  of  the  0(  P)  +  H20(  A-j )  surfaces. 

Assuming  a  self-consistent  field  (SCF)  starting  point,  the 
defining  equations  of  MBPT  are 


HQ  =  l  [h(i)  +  u ( i ) ] 


n  ,  n 

v  5  J.  rij  - 1  u<*> 

i'J  i 


Hj  -  E 
0*0  0  0 


Eo  =  I  £f 


t>0  =  A  (xi(1)x2(2).  .  -xn(n)) 


E  =  E0  *  7  ♦0I»[(E0-H0)-,V1PU  >L 
P=  ° 
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The  unperturbed  hamil Ionian,  HQ,  consists  of  the  usual  one-electron 
hamiltonian  and  the  SCF  effective  one-particle  potential.  4Q  is  the  SCF 
function  composed  of  the  occupied  orbitals  (x^l  associated  with  the 
orbital  energies  {e^l.  The  linked  diagram  theorem,  Eg.  (6),  defines  the 
energy  as  a  sum  of  diagrams  where  the  subscript  L  restricts  them  to  "linked" 
diagrams.  Each  diagram  provides  an  algebraic  expression  in  terms  of  the 

transformed  two-electron  integrals,  <xpxq ! I xrxs>  =  <pq!|rs>  =  (  xp>  rlxqxs) 

-  ( x  Xs i XqXr ) »  and  the  orbital  energies  le.i  .  By  simply  evaluating  these 
expressions  composed  solely  of  information  generated  in  the  SCF  calculations, 
the  correlation  contributions,  En,  for  n_>2  are  obtained.  This  is  a  very 
efficient  process  leading  to  easily  vectorizable  codes  that  are  ideally 
suited  to  the  new  generation  of  computers. 

As  an  example,  in  second-order 

E2  •  X  !<1JI|,bM2X  i  *  '  'a  "  ‘b  1  <7) 

a  b 

where  i,  j  are  orbitals  occupied  in  4^  and  a,  b  are  excited  SCF  orbitals. 
This  term,  which  takes  0  negligible  amount  of  time  in  addition  to  the 
SCF  calculation,  accounts  for  the  bulk  of  the  correlation  effects  in 
a  problem.  However,  we  have  found  that  to  get  very  good  relative  enerqy 
differences  at  adjacent  points  it  is  usually  necessary  to  include  E^  and 
most  of  E4.(6,7)  The  latter  model,  SDQ-MBPT(4).  refers  to  oil  single,  double, 
and  quadruple  excitation  effects  through  fourth-order.  This  model  does 
neglect  the  triple  excitations,  but  their  inclusion  would  excessively 
increase  the  required  computer  time  for  what  appears  to  be  little 


additional  benefit.  In  comparisons  with  the  infinite-order  coupled- 
cluster  doubles  model,  one  seldom  finds  more  than  1  kcal/mole  difference 
due  to  the  hiqher  than  fourth-order  contributions,^  so  this  model  is 
highly  stable  with  respect  to  possible  higher-order  corrections  of  the 
same  type  as  those  included  in  fourth-order!^  1 

A  number  of  applications  have  shown  that  calculations  at  the 
SDQ-MBPT(4)  level  and  with  at  least  a  double-zeta  plus  polarization 
basis  set  are  capable  c'  predicting  the  binding  energies  in  the  reactions 
ZBH^B^Hg,  BH3  +  CCMH^BCO,  and  BH^  +  NHj-^BNH^  to  within  +2  kcal/mole.^ 
Similarly,  the  activation  barriers  and  tne  isomerization  energies  of  the 
unimolecular  isomerization  reactions  CH^NC-^CH^CN  and  HNC-HCN^  are 
exceptionally  accurate.  In  applications  to  the  potential  energy  surface 
for  the  reaction  HCO-H+CO; 1 ^  the  quartic  force  field  of  ^o!8)  and 
the  step-wise  decomposition  of  methanol !^  it  is  shown  that  methods 
at  this  level  are  reliable  for  determining  bond  lengths  and 
bond  angles,  and  in  the  case  of  i^O,  for  calculating  force-constants 
even  including  the  higher-order  anharmonic  constants.  One  therefore  has 
every  reason  to  expect  that  the  SDQ-MBPT(4)  model  with  an  adequate  basis 
set  should  provide  an  accurate  description  of  the  ground  state 
triplet  surface. 
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Ill,  0(3P)  ♦  H20(1A1 )  COLLISION 
A.  Potential  Energy  Surface 

0  1 

Before  considering  the  global  aspects  of  the  0(  P)  +  H^0(  A-j ; 
surface,  it  is  pertinent  to  focus  on  the  two  individual  fragments.  The 
SCF  and  correlation  energies  for  a  range  of  basis  sets  are  presented  in 
Table  I. 

The  (5s4p2d/3slp)  Slater  orbital  (STO)  basis  has  been  used 

(o) 

in  a  study  of  the  quartic  force  field  of  H^O.  These  calculations 
are  done  using  SCF,  SD-CI,  and  with  various  many-body  models.  The 
theoretical  results  agree  well  with  the  experimental  force  constants, 
with  the  many-body  values  showing  better  agreement  than  SD-CI.  In  fact, 
the  most  recent  experimental  and  normal  coordinate  study  revises  several 

of  the  previously  accepted  force  constants  to  be  more  consistent  with 

(81 

these  theoretical  calculations.  All  of  the  cubic  and  quartic  constants 
are  determined  from  the  calculations,  while  the  experimental  studies 
typically  set  many  of  these  constants  to  zero  to  facilitate  obtaining 
the  force  field  from  the  normal  coordinate  analysis.  Hence,  the 
computed  SDQ-MBPT ( 4 )  quartic  force  field  for  H-,0  is  used  as  input  to 
the  local  description  of  H90  in  the  0  +  H90  surface  fit.  This  is  described 
in  detail  in  Section  III.B. 

Since  it  would  be  very  expensive  to  use  a  STO  basis  far  the 
four  atom  surface  calculations,  a  double-zeta  plus  polarization  (DZP) 
contracted  Gaussian  basis  is  used  to  determine  the  behavior  of  the  surface 
describing  the  approach  of  the  oxygen  atom  to  the  H^O  molecule.  The 
results  in  this  basis  are  compared  fo  H^O  to  the  STO  basis  and  other 
contracted  Gaussian  bases'  in  Table  I.  In  the  case  of  the  DZP  results. 
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TABLE  1.  ENERGIES  FOR  0(3P)  AND  H^Of  > At ) ,  AND  GEOMETRIES  FOR  HpO(^Ai) 
IN  DIFFERENT  BASI$  SETS.  ALL  ENERGIES  ARE  IN  HARTPEE  A.lj.  P 
DISTANCES  IN  BOORS 


The  double-zeta  plus  polarization  basis  are  those  used  in  the  complete 
surface  study  of  the  0(3p)  +  H^OpAi)  surface.  All  calculations  are  done 
by  freezing  the  Is2  core  in  the  oxygen  atoms  at  the  SCF  level.  Hence, 
the  correlation  energies  reported  in  this  case  are  solely  for  the  valence 
electrons. 

b  P.  S.  Bagus,  B.  Liu,  and  H.  F.  Schaefer,  III,  Phys.  Rev.  A2,  555  (1970). 

c  R.  J.  Bartlett,  I.  Shavitt,  and  G.  D.  Purvis,  J.  Chem.  Phys.  71_,  281  (1979). 

^  Valence  shell  correlation  energies  are  obtained  from  the  experimental 
values  by  subtracting  the  KK,  and  KL  shell  contributions  to  the 
correlation  energy  obtained  by  R.  K.  Nesbet,  Phys.  Rev.  U5,  2  (1961), 
for  the  He  isoelectronic  sequence. 

e  J.  A.  Pople  and  J.  S.  Binkley,  Mol.  Phys.  29,  599  (1975). 

f  A.  R.  Hoy  and  P.  R.  Bunker,  J.  Mol.  Spectrosc.  74,  1  (1979). 

9  A.  R.  Hoy,  I.  M.  Mills,  and  G.  Strey,  Mol.  Phys.  24,  1265  (1972). 
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the  oxygen  Is  electrons  are  kept  frozen  at  the  SCF  level,  so  only  the 
valence  correlation  energy  is  computed.  The  other  results  are  modified 
to  provide  a  valence  correlation  result  for  comparison. 

Although  the  best  calculation  provides  about  80%  of  the 
valence  correlation  energy  the  DZP  basis  still  gives  69%.  In  addition, 
the  bond  length  and  bond  angle  in  this  basis  differ  only  slightly  from 
the  STO  calculation.  These  results  suggest  that  we  should  adequately 
describe  the  shape  of  the  surface  at  this  level.  In  fact,  as  will  be 
seen  below,  the  0  +  Ho0  triplet  surfaces  are  sufficiently  simple,  that 

4_ 

even  a  DZ  basis  would  not  lead  to  much  of  an  error  in  the  surface.  For 
the  oxygen  atom,  the  DZP  basis  accounts  for  61"  of  the  valence  correlation 
energy,  leading  to  a  separated  limit  of  -151.17567  a.u.  One  additional 
piece  of  information  provided  in  Table  1  is  tne  barrier  to  inversion  for 
H^O.  This  is  found  to  be  1.5eV  by  optimizing  the  nonlinear  and  linear 
HoO  structures.  Since  for  coll'sion  energies  of  0.5  to  2.5  eV  inversion  can 
occur,  the  cross-sections  will  be  sensitive  to  the  possibility  of  inver¬ 
sion  in  the  bending  mode. 

The  collision  between  0{  p)  atoms  and  ground  state  H^O  results 
in  three  electronic  surfaces  that  correlate  with  0(^P).  These  are 
illustrated  in  the  correlation  diagram  in  Figure  1,  for  different 
directions  of  approach  of  the  oxygen  atom.  The  one  pertinent  reaction 
channel  giving  20H(2i:)  molecules  is  also  indicated.  The  p-imary 
objective  of  the  present  work  is  to  investigate  the  lowest  of  these 
surfaces.  Since  all  three  surfaces  are  relevant  to  a  description  of 
the  collision  dynamics,  the  higher  two  surfaces  will  be  considered  to 
some  degree. 
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Depending  upon  the  symmetry  involved  when  0(  P)  approaches  H^O, 
the  three  electronic  states  may  all  be  of  different  symmetries,  two  different 
symmetries,  or  all  of  the  same  symmetry.  As  long  as  the  symmetries  are 
different,  the  UHF-MBPT  approach  can  be  used  to  obtain  the  surface. 

However,  it  cannot  conveniently  calculate  an  excited  surface  of  the  same 
symmetry  as  a  lower  state,  since  the  UHF  solutic  is  difficult  to  obtain 
and  may  suffer  from  large  amounts  of  spin  contamination. 

If  H20  is  placed  in  the  yz  plane,  such  that  the  z-axis  corresponds 

to  the  two-fold  rotation  axis,  as  in  Figure  2,  different  approaches  may  be 

considered.  An  approach  of  the  colliding  atom  in  the  yz  plane  along  the 

molecular  axis  maintains  the  C2v  syrmetry  of  the  system.  For  this  approach 

13  3 

it  is  possible  to  obtain  B-j ,  B^,  and  A2  surfaces.  Similarly,  an 

approach  maintaining  a  single  plane  of  symmetry,  such  as  along  the  y-axis 
but  in  the  plane  of  the  molecule,  results  in  C$  symmetry  with  one  A* 
and  two  A"  states.  If  there  is  no  plane  of  symmetry,  only  the  single 
lowest  state  is  obtained.  Figure  2  shows  the  distribution  of  ab  initio 
calculations  for  H20  at  equilibrium.  Also  calculations  using  various 
displacements  of  H20  are  used.  To  date,  more  than  90  points  on  the 
various  surfaces  have  been  computed.  Some  of  these  are  listed  in  Tables  2-9. 

Figures  3-7  show  curves  describing  different  cuts  through  the 
triplet  surfaces.  It  is  apparent  from  all  of  the  figures  that  these  are 
essentially  repulsive,  at  least  for  small  0-H  bond  displacements.  This 
should  be  contrasted  with  a  surface  describing  the  interaction  between 
H20  and  a  singlet  oxygen  atom  which  would  have  a  number  of  local  extrema 
including  the  three  minima  in  the  H207  (hydrogen  proxide)  structure. 
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TABLE  h20  AT  EQUILIBRIUM.3  COLLIDING  O-ATOM 
APPROACHING  ALONG  THE  NEGATIVE  Z-AXIS 


Distance  to  0  Atom(bohr)^  SC F  SDQ-MBPT(4) 


^B-i  Ground  State 

-2.6 

-150.56647 

-150.93642 

O 

m 

1 

-150.72168 

-151.06774 

-4.0 

-150.83686 

-151.16515 

-6.0 

-150.85317 

-151.17711 

-8.0 

-150.85259 

-151.17621 

^  First  Excited  State 

-2.5 

-150.52779 

-150.88583 

-3.0 

-150.70554 

-151.04874 

-4.0 

-150.83389 

-151.16192 

-6.0 

-150.85288 

-151.17682 

-8.0 

-15C. 85253 

-151.17615 

A2  Second  Excited  State 

-2.5 

-150.28980 

-150.64314 

-3.0 

-150.57280 

-150.91766 

-4.0 

-150.80587 

-151.1  3298 

-6.0 

-150.84Q93 

-151.17386 

t 

O 

00 

1 

-150.85171 

-151.17530 

a  Coordinates  of  H^O  molecule  in  Dohrs  are  (  C  :0.0,0.0,-P  J 2371  ;H:O.P  ,+l  .43035, 
0.98168).  Total  energies  in  Hartrees. 
b  Distance  is  measured  from  the  center-of-mass  of  the  H^O  molecule. 
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TABLE  3  .  H?0  AT  EQUILIBRIUM.3  COLLIDING  O-ATOM 

APPROACHING  ALONG  POSITIVE  Z-AXIS 


Distance  to  0  Atom  (bohr)b  SCF  SDQ-MBPT(4) 


3B^  Ground  State 

+2.5 

-150.47027 

-150.79727 

+3.0 

-150.69253 

-151.02192 

+4.0 

-150.82630 

-151.15808 

+6.0 

-150.85192 

-151.17593 

+8.0 

-150.85212 

-151.17578 

3B2  First  Excited  State 


+2.5 

-150.43655 

-150.79260 

+3.0 

-150.67233 

-151.02:65 

+3.5 

-150.77522 

-151.11363 

+4.0 

-150.82018 

-151.15090 

+4.5 

-150.83914 

-151.16630 

+6.0 

-150.85111 

-151.17532 

+8.0 

-150.85196 

-151.17562 

3A2  Second  Excited  State 

+2.5 

-150.20481 

-150.58459 

+3.0 

-150.56853 

-150.91886 

+4.0 

-150.80486 

-151.13336 

+6.0 

-150.85295 

-151.17704 

+8.0 

-150.85282 

-151.17640 

a  Coordinates  of  H20  molecule  in  bohrs  are( 
0.98161)-  Total  energies  in  Hartrees. 
b  Distance  is  measured  from  center-of-mass 


0:0.0,0.0,-0.12371 ;  H:0.0, ±143085, 
of  the  H^O  molecule. 
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TABLE  4 


H20  at  equilibrium-3  colliding  o-atom  approaching 

ALONG  Y-AXIS 


Distance  to  O-Atom(bohr) 

SCF 

SDQ-MBPT(4) 

Ground  State 

+2.5 

-150.38065 

-150.75549 

+3.0 

-150.65074 

-151.00760 

+3.5 

•  150.76810 

-151.10805 

+4.0 

-150.82002 

-151.15058 

+5.0 

-150.85025 

-151.17531 

+6.0 

-150.85376 

-151.17801 

+8.0 

-150.85215 

-151.17829 

A"  First  Excited  State 

+2.5 

-150.33535 

-150.70048 

+3.0 

-150.63058 

-150.98266 

+3.5 

-150.75750 

-151 .09635 

+4.0 

-150.81343 

-151.14387 

+4.5 

-150.83676 

-151.16344 

+5.0 

-150.84609 

-151 .17116 

+6.0 

-150.85112 

-151.17538 

+8.0 

-150.85209 

-151.17825 

a  Coordinates  of  H?0  molecule  in  bohrs  are  (  0*0.0,0.0,-0.12371 ; 

H:0.0,+1 .43085,0.98168). 

Total  energies  in 

Hartrees . 

b  Distance  is  measured  from  the  center-of-mass  of  the  HgO  molecule. 
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TABLE  5  .  H20  AT  EQUILIBRIUM.9  COLLIDING  O-ATOM 

APPROACHING  ALONG  THE  X-AXIS 


Distance  to 

j-ator?/',ohr)k  SCF 

SDQ-MBPT(4) 

A'  Ground  State  (Approach  in  direction  0°) 

3.0 

-150.72692 

-151 .07130 

4.0 

-150.83585 

-151.16543 

8.0 

-150.8523/ 

-151.17565 

A'  Ground  State  (Approach  in  direction  +45°) 

3.0 

-150.70467 

4.0 

-150.82954 

-151.15960 

6.0 

-150.85286 

-151 .17663 

A'  Ground  State  (Approach  in 

direction  -45°) 

3.0 

-150.71522 

-151.05884 

4.0 

-150.83500 

-151 .16349 

6.0 

-150.85003 

-151 .17373 

Coordinates  of  H20  molecule  in  bohrs  are  (  0:0.0,0.0,-0.12371; 

H: 0.0, +1 .43085,  0.98164).  Total  energies  in  Hartrees. 

Distance  is  measured  from  the  center-of-mass  of  the  H20  molecules. 


27 


TABLE  6.  H?0  IN  A  LINEAR  CONFIGURATION.3 

COLLIDING  O-ATOM  APPROACHING  ALONG  THE  Z-AXIS 


Distance  to  0- 

-atom  (bohr)*3  SCF 

SDQ-MBPT(4) 

3 

B.|  Ground  State 

-2.5 

-150.5223 

-150.8949 

-3.0 

-150.6875 

-151.0306 

-4.0 

-510.7869 

-151.1127 

-6.0 

-150.8017 

-151.1216 

h  First  Excited  State 

-2.5 

-150.4982 

-150.8603 

-3.0 

-150.6720 

-150.0123 

-4.0 

-150.7826 

-151.1079 

-6.0 

-150.8011 

-151.1211 

3 

A^  Second  Excited  State 

-2.5 

-150.3062 

-150.6480 

-3.0 

-150.5444 

- 

-4.0 

-150.7556 

-151.0779 

-6.0 

-150.7991 

-151.1188 

r--  -  t.  i  ■  sr  ,i.«i  X  r... i^.1  1 . 

a  Coordinates  of  HjO  in  bohrs  are  (0  :  0.0, 0.0, 0,0;  H:0.0,  +1.735236,  0.0). 
Total  energies  in  Hartree  atomic  units. 

h  Distance  is  measured  from  the  0-atom  in  the  linear  H^O  configuration. 
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TABLE  7.  H20  IN  ASYMMETRIC  STRETCH.3  COLLIDING 

O-ATOM  APPROACHING  ALONG  THE  Y-AXIS 


Distance  to  0-Atom  (bohr)^ 

SCR 

SDQ-MBPT(4) 

A'  Ground  State 

+4.0 

-150.81538 

-151.14604 

+3.5 

-150,76315 

-151.10243 

+3.0 

-150.64065 

-150.99692 

+2.5 

-150.35848 

-150.73519 

-2.5 

-1 50.37837 

-150.75(94 

-3.0 

-150.64653 

-151.00527 

-3.5 

-150.76036 

-151.10242 

-4.0 

-150.81242 

-151.14442 

-6.0 

-150.84798 

-151.17208 

-8.0 

-150.84703 

-151.17050 

a  Coordinates  of  H^O  molecule 

in  bohrs  are  (0 

:  0.0,0.0,-0.12371, 

HI:  0.0,1.3640,  0.90480,  H2: 

:  0.0,  -1.52390 

,  1.02520). 

Total  energies  in  Hartrees  atomic  units. 

Distance  is  measured  from  the  center-of-mass  of  the  H^O  molecule  at 
its  original  equilibrium  position. 
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TABLE  8  .  H?0  IN  A  BENT  CONFIGURATION.8  COLLIDING 

O-ATOM  IS  APPROACHING  ALONG  THE  Z-AXIS 


Distance  to  0-Atom(bohr)b  SCF  SDQ-MBPT(4) 


3 

B.|  Ground  State 


-2.5 

-150.57023 

-150.95284 

-:.o 

-150.72326 

-151.06877 

-3.5 

-150.80442 

-151.13847 

+2.5 

-150.48542 

-150.85942 

+3.0 

-150.70036 

- 

+3.5 

-150.78971 

-151.12909 

Coordinates  of  H^O  molecule  in  bohrs  are  (  0:0.0,0.0,-0.10968; 

HI : 0.0+1 .52650,  0.87032).  Total  energies  in  Hartrees  atomic  units. 
Distance  to  the  0-atom  is  measured  from  the  center-of-mass  of  H^O  i 
original  equilibrium  configuration. 
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Since  the  ground  state  surface  is  largely  repulsive  the  number  of  geometries 
sampled  in  these  calculations  is  not  as  large  as  would  be  necessary  for 
a  surface  with  local  extrema. 

In  Figure  3  are  shown  the  three  surfaces  obtained  at  the  SCF 

(i.e.,  UHF)  level.  The  strong  repulsive  character  is  evident.  Also, 

3  3 

there  is  only  a  modest  difference  in  the  and  curves,  with  slightly 
more  distinction  from  the  A£  curve.  As  shown  in  Figure  4,  the  behavior 
is  unchanged  at  the  SDQ-MBPT(4)  level.  In  the  case  of  the  correlated 
treatment,  the  total  energy  is  found  to  be  much  lower  and  the  curves  are 
displaced  accordingly,  but  from  the  viewpoint  of  collisir  ,  dynamics, 
cross  sections  will  be  affected  only  by  the  shape  of  the  curves.  A  closer 
comparison  demonstrates  that  the  SC'r  and  SDQ-MBPT(4)  curves  for  H^O 
at  equilibrium  are  almost  superimposable. 

This  behavior  is  further  emphasized  in  Figure  5  when  the  SDO- 
MBPT(4)  and  SCF  curves  are  compared  for  an  approach  along  the  y-axis  of 
the  H^O  molecule.  Although  the  difference  between  the  correlated  and 
uncorrelated  curves  is  not  constant  for  this  approach,  it  is  close  to  being 
so.  The  correlation  energy  is  actually  more  nearly  constant  for  collisions 
approaching  along  the  z-axis. 

In  Figure  6  comparisons  of  the  ^  state  at  four  different  levels 
of  approximation  are  presented.  In  addition  to  the  similarity  in  shape  of 
the  SCF  and  SDQ-MBPT(4)  curves,  it  is  also  evident  that  a  restriction  of 
the  basis  set  to  DZ  (double-zeta )  instead  of  tne  polarized  DZP 
basis  also  has  essentially  no  effect  on  the  shape.  A  superposition  of 
these  curves  would  show  no  significant  differences. 
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A  State 


-150.40 


Bohrs 


FIGURE  6.  FOUR  DIFFERENT  APPROXIMATIONS  TO  THE  3B?  STATE 
FOR  0(3P)  COLLIDING  WITH  H?0  AT  EQUILIBRIUM 
APPROACHING  ALONG  THE  POSITIVE  Z-AXIS.  DIFFERENT 
BASIS  SETS,  DOUBLE-ZETA(DZ)  AND  DOUBLE-ZETA  PLUS 
POLARIZATION  (DZP),  AT  BOTH  THE  SCF  AND  SDQ-MBPT(41 
LEVEL  ARE  ILLUSTRATED. 
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This  observation  suggests  two  important  points.  The  theoretical 
calculations  at  several  different  degrees  of  sophistication  show  very 
little  differences,  so  that  convergence  to  the  correct  potential  surface 
has  been  effectively  achieved.  Further  additions  to  the  basis  set  or 
to  the  correlation  effects  included  would  not  appear  to  cause  important 
changes  in  the  surface.  Also,  additional  work  at  the  m  ch  less  expensive 
double-zeta  level  should  therefore  provide  useful  results.  This  is  an 
important  observation  that  permits  using  a  DZ  basis  for  the  more  expensive 
0(3P)  +  C02(]Zg)  PES. 

3 

Another  benefit  of  the  observed  repulsive  nature  of  the  0(  P)  +  H^O 
surface  is  that  the  analytical  representation  of  the  surface  is  much  simpler 
than  when  minima  must  be  reproduced. 

Figures  7  and  8  illustrate  the  shape  of  the  surface  for  two 
possible  displacements  of  the  H^O  molecule  from  equilibrium,  the  linear 
configuration  and  an  asymmetric  stretch.  Again  the  overall  shapes  of 
the  curves  are  similar,  but  there  are  observable  differences  caused  by 
the  displacement  of  b^O  as  a  function  of  0-atom  position  which  must  be  realist 
cally  described  in  order  for  the  dynamical  calculations  to  result  in 
reliable  cross  sections. 

Probably  the  distorted  structures  having  the  most  significant 
effect  on  the  001  cross  sections  involve  asymmetric  stretch  displacements. 

3 

This  is  also  expected  to  be  the  route  taken  for  the  reaction  0(  P)  + 

1  2 

H20(  A-] )  •>  20H(  n),  as  shown  in  Figure  1.  Since  the  bond  energy  of 
H  -  OH  is  119  kcal/mole  and  0-H  is  102.4,  this  reaction  is  endoergic  by 
only  17  kcal/mole.  At  the  collision  velocities  of  interest  in  the  current 
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work,  it  does  not  appear  likely  that  the  reaction  itself  will  withdraw 
much  flux  from  the  vibrational  excitation  channels,  since  little  time 
is  available  to  sample  the  reactive  channel.  However,  for  those 
trajectories  in  which  large  displacements  of  the  OH  bond  occur,  the 
potential  surface  used  in  this  work  may  not  be  entirely  accurate.  Even 
though  a  number  of  points  involving  asymmetric  stretches  of  H^O  have 
been  computed,  it  is  outside  the  scope  of  the  present  work  to  consider 
reactive  collisions.  Hence,  there  are  not  enough  points  computed 
to  locate  and  classify  a  possible  transition  state  region.  This  should 
be  done  in  future  work. 
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B.  The  Analytical  Representation  of  the 
0 C3P )  +  H20  Potential  Energy' Surface 

The  accurate  representation  of  multidimensional  potential 
energy  surfaces  has  always  been  a  formidable  problem.  A  common 
approach  is  to  employ  an  expression  that  incorporates  as  much  physical 
insight  as  possible  in  the  functional  form,  and  which  has  a  number  of 
parameters  that  are  adjusted  to  empirical  data.  Examples  of  this 
approach  are  found  in  applications  of  the  London-Eyring-Polyani-Sato 
(LEPS),  valence-bond  (VB),  diatomics-in-molecules  (DIM),  and  many-body 
expansion  methods  for  polyatomic  systems!^'^8^ 

Recent  developments  in  ab  initio  methods  for  calculating 
potential  energy  surfaces  will  require  new,  systematic  approaches  for 
obtaining  useful  analytical  representations.  Progress  has  been  made  in 
this  area  using, for  example,  mul tidimensional  splines^’^,  polynominal 

(21 )  (2?  23) 

root  interpolation,  and  many-body  expansions.  ’  An  important  goal 

is  to  develop  methods  that  can  be  systematically  applied  and  that  require 

a  relatively  small  number  of  ab  initio  points 

The  method  used  in  the  present  work  is  based  on  the  Sorbi e-Murrell 

(24 ) 

many-body  representation  of  a  potential  energy  surface,  '  which  has  been 

(18) 

applied  by  Murrell  and  coworkers  to  a  number  of  systems.  1  For  studying 
non-reactive  collisional  excitation  processes  involving  an  oxygen  atom 
colliding  with  a  water  molecule,  the  many-body  expansion  can  be  written 

V(0,  h20)  =  v00(RAB)  +  v0H(RAC)  +  v0H  (RAD) 

+  vh2o  (RBC’  RBD’  W  +  V0HH^RAC’  rad* 
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+  vooh^rab’  rbc’  rac^  +  vooh^rab’  rbd*  rad} 


+  V0(WRAB’  rac*  rad*  rbc*  rbd*  rcd^ 


In  this  expression  A  and  B  represent  the  reagent  atom  and  oxygen  atom  in 

the  molecule,  respectively,  while  C  and  D  represent  che  hydrogen  atoms. 

Since  reactive  channels  are  excluded  from  the  present  study, 

one  can  choose  the  fragment  potentials  to  ensure  that  the  contribution 

of  the  four-body  term  is  small.  Since  the  ab  initio  points  show  a 

repulsive  surface  (in  the  region  sampled),  the  two-body  terms  are 

chosen  to  be  exponentials.  The  term  g  is  taken  to  be  an  accurate 

^  ( 8  ?5 1 

ground-state  surface  for  the  water  molecule.  *  'All  of  the  remaining 
three-  and  four-body  terms  are  represented  as  products  of  polynominals 
times  suitable  cutoff  functions. 


The  choice  of  exponentials  for  the  intermolecular  pair  potentials 
is  made  so  that  the  simple  sum-of-pairs  plus  intramolecular  force  field 
representation 


v^(o,h2o)  =  vqo^rab^  +  voh  ^rac^  +  VW  +  vh2o^rbc,rbd,rcd^ 

contains  the  essential  features  of  the  surface.  The  two-body  terms  are 

represented  as  sums  of  two  tAponential  functions 

-a,R  -a7R  r 

vab(R)  ■  A]  e  +  A2  e  -C/Rb 


with  the  exponents  adjusted  to  the  slope  of  the  surface  along  two  cuts. 

The  linear  coefficients  are  determined  by  linear  regression  using  a 
large  number  of  points  for  which  H20  is  at  its  equilibrium.  The  resulting 
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fit  is  tnen  corrected  by  adding  in  the  additional  three-body  and  then 
four-body  terms,  using  linear  regression  at  each  stage  to  determine  the 
polynomial  coefficients.  This  time,  points  that  represent  vibrational 
displacements  of  h^O  away  from  equilibrium  art  included.  This  allows 
the  fit  to  properly  describe  the  variation  of  the  force  constants  of  the 
water  molecule  with  the  approach  of  the  oxygen  atom.  It  should  be 
mentioned  that  sequential  inclusion  of  the  higher  order  terms  is  necessary 
to  avoid  unphysical  behavior  of  the  fit  in  regions  where  ab  initio  points 
are  not  available.  This  undesirable  behavior  always  occurs  if  all  of 
the  coefficients  are  adjusted  simultaneously. 

The  decay  functions  used  to  damp  out  the  polynomials  are  of 
the  form^  ^ 

N 

II  [1  -  tanh  ( Y i S . /2 ) 3  , 
i  =  l  1 

where  N  is  3  or  6,  and  =  R ^ - R ^ Q .  Both  R.  and  y,.  are  chosen  to  . 
the  functions  so  that  the  ab  initio  points  could  be  fit  without  the  polynomial 
coefficients  becoming  too  large.  This  helps  in  smoothing  the  overall 
representation  of  the  surface. 

A  very  important  feature  of  the  many-body  approach  is  the 
convenience  with  which  information  about  the  fragments  can  be  included. 

Since  this  study  emphasizes  collisional  energy  transfer,  it  is  important 
to  have  a  reliable  surface  for  the  water  molecule  that  reproduces  the 
known  force  constants. 

The  first  potential  used  was  a  very  accurate  quartic  force 
field  obtained  by  SDQ-MBPT(4) . v  ’  Because  high  velocity  collisions 
are  found  to  greatly  perturb  the  bending  mode,  a  potential  that  is 


reliable  further  from  equilibrium  is  required.  Consequently,  the  Sorbie- 
Murrell  representation  of  the  water  molecule  is  adopted,  but  with  parameters 

(25) 

determined  so  that  the  fit  reproduced  the  MBPT  quartic  force  field 
rather  than  experiment. 

In  Table  10  a  comparison  is  made  of  the  force  constants  obtained 

(25) 

in  Reference  8  with  those  of  the  fit,'  and  the  agreement  is  seen  to  be 
excellent.  As  seen  in  Figure  9,  however,  the  bending  potential  from  the 
fit  for  Ho0  is  much  more  realistic  away  from  equilibrium  than  is  the 
quartic  force  field.  This  is  demonstrated  by  the  existence  of  a  saddle 
point  for  the  bending  motion  for  a  linear  geometry  of  h^O  for  the  many- 
body  fit  with  a  barrier  of  1.296  eV  above  the  water  equilibrium  energy, 
while  the  quartic  force  field  does  not  even  have  a  saddle  point. 


Experimental  estimates  are  that  the  barrier  is  1.37  eV  above  the  equilibrium 

(19) 

energy,'  while  calculations  using  the  present  DZP  basis  give  1.5  eV. 

Using  the  many-body  fit  for  results  in  an  order  of  magnitude  increase 
in  the  asymmetric  stretch  excitation  cross  section  at  1.5  eV  relative 
energy.  It  would  appear  then  that  the  choice  of  V^q  is  critical  for 
studies  of  collisional  energy  transfer  in  water.  Contours  for  the 
stretching  displacements  of  H^O  are  in  Figure  10. 

The  ab  initio  points  are  selected  with  the  dynamics  applica¬ 
tions  in  mind.  Since  high-velocity  collisional  excitation  of  the 


vibrational  modes  is  of  interest,  most  of  the  calculations  are  made 
for  0(  P)  in  the  plane  of  the  E^O  molecule.  In  Figure  2  the  distribution 
of  points  in  this  plane  is  shown  for  E^O  at  its  equilibrium  qeometry. 

As  discussed  previously  calculations  are  also  made  for  cut-of-plane 


geometries  and  for  nonequilibrium  E^O  configurations. 
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reliable  further  from  equilibrium  Is  required.  Consequently,  the  Sorbie- 
Murrell  representation  of  the  water  molecule  is  adopted,  but  with  parameters 

determined  so  that  the  fit  reproduced  the  MBPT  quartic  force  field' 

(24) 

rather  than  experiment. 


In  Table  10  a  comparison  is  made  of  the  force  constants  obtained 

(25) 

in  Reference  8  with  those  of  the  fit,  and  the  agreement  is  seen  to  be 

excellent.  As  seen  in  Figure  9,  however,  the  bending  potential  from  the 

fit  for  H^O  is  much  more  realistic  away  from  equilibrium  than  is  the 

quartic  force  field.  This  is  demonstrated  by  the  existence  of  a  saddle 

point  for  the  bending  motion  for  a  linear  geometry  of  h^O  for  the  many- 

body  fit  with  a  barrier  of  1.296  eV  above  the  water  equilibrium  energy, 

while  the  quartic  force  field  does  not  even  have  a  saddle  point. 

Experimental  estimates  are  that  the  barrier  is  1.37  eV  above  the  equilibrium 
(19) 

energy,  while  calculations  using  the  present  DZP  basis  give  1.5  eV. 

Usinq  the  many-body  fit  for  V^q  results  in  an  order  of  magnitude  increase 
in  the  asymmetric  stretch  excitation  cross  section  at  1.5  eV  relative 
energy.  It  would  appear  then  that  the  choice  of  VH  q  is  critical  for 
studies  of  collisional  energy  transfer  in  water.  Contours  for  the 
stretching  displacements  of  H20  are  in  Figure  10. 

The  ab  ini ti o  points  are  selected  with  the  dynamics  applica¬ 
tions  in  mind.  Since  high-velocity  collisional  excitation  of  the 
vibrational  modes  is  of  interest,  most  of  the  calculations  are  made 
for  O('T)  in  the  plane  of  the  H^O  molecule.  In  Figure  2  the  distribution 
of  points  in  this  plane  is  shown  for  H20  at  its  equilibrium  geometry. 

As  discussed  previously  calculations  are  also  made  for  out-of-plane 
geometries  and  for  nonequilibrium  H^O  configurations. 
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TABLE  10.  FORCE  CONSTANTS  FOR  H?0 


EXPa 

MBPTb 

FITC 

u  1 

3831.5 

3865.0 

3864.9 

U)^ 

1648.8 

1687.4 

1687.4 

u>3 

3942.2 

3975.0 

3975.0 

km 

-302.5 

-304.2 

-304.2 

k222 

63.6 

42.9 

42.9 

kl  22 

167.4 

148.6 

148.6 

k21 1 

-53.1 

-61.8 

-61 .8 

kl  33 

-927.8 

-914.1 

-914.1 

k233 

-138.8 

-111.7 

-111.6 

'llll 

31.9 

31.5 

31.5 

2222 

2.1 

-2.6 

-2.6 

'3333 

35.4 

32.0 

32.0 

'1122 

-85.6 

-75.1 

-74.8 

'1133 

201.3 

190.3 

190.4 

k2233 

-101.1 

-91.7 

-91 .4 

+  - 1 
Units  are  cm  . 

d  Hoy-Mi  11 s-Strey. 
b  MBPT 

c  Many-Body  Fit. 
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Potential  Energy  (eV) 


0-H  Bond  Distance,  A 


FIGURE  10.  POTENTIAL  CONTOURS  FOR  STRETCHING  DISPLACEMENTS 
FOR  Ho0.  CONTOUR  INCREMENT  IS  1  eV 


A  number  of  fits  to  the  computed  ab  initio  points  are  made, 
using  various  combinations  of  points  and  terms  in  the  many-body  expansion. 
Those  that  are  discussed  here  are  used  in  trajectory  calculations  that 
will  be  presented  in  a  later  section.  The  various  fits  are  denoted  by 
combinations  of  alphanumeric  characters,  such  as  E2Q.  This  fit  is 
determined  using  only  two-body  terms  for  the  intermolecular  potential. 

uses  three-body  terms  to  second  order,  and  F^g  uses  four-body  terms 
to  second  order.  E^^uses  three-body  terms  to  fourth  order.  E^  ,  F^, 
and  differ  essentially  in  the  weights  used  for  the  linear  regression, 
and  the  choice  of  parameters  for  the  cutoff  functions. 

In  Figure  Hare  shown  bending  potentials  for  motion  of  the 
incominq  0-atom  in  the  plane  of  the  1^0  molecule  at  a  constant  distance 
of  3.0  flQ  from  the  center-of-mass  of  h^O.  0r  corresponds  in  Figures  11-13 
to  an  approach  along  the  z-axis  toward  the  two  hydrogen  atoms,  and  180° 
corresponds  to  an  approach  toward  the  0-atom  end  of  the  molecule.  The 
two  curves  labeled  F  represent  the  best  fits,  although  the  two-body 
fit  E2a  is  realistic.  Subtle  changes  in  the  choice  of  weights  or  of 
the  nonlinear  parameters  can  significantly  affect  the  quality  of  the  fit. 
The  three-body  fits  in  Figure  11  involve  only  second-order  polynomial  terms. 

The  influence  of  higher-order  terms  is  illustrated  in  Figure  12  . 
Fit  E.j^  employs  three-body  terms  to  fourth  order,  and  provides  a  fit 
to  most  of  the  74  points  used  to  within  aoout  1%.  Although  this  is  an 
accurate  fit  to  the  computed  points,  the  possibility  of  the  occurrence 
of  unphysical  oscillations  prevents  its  use  in  trajectory  studies  as  many 
of  the  polynomial  coefficients  are  quite  large.  Further  work  at  this 
leve1  would  be  likely  to  produce  a  useful  surface. 
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FIGURE  II.  BENDING  PO 7EN71AL  FOR  0(3P)  ABOUT  HoO  IN  THE  PLANF 
OF  THE  MOLECULE  AT  A  CONSTANT  DISTANCE  OF  3.5  aQ 
FROM  THE  MOLECULAR  CENTER-OF-MASS.  THE  VARIOUS 
FITS  ARE  DISCUSSED  IN  THE  TEXT.  COMPARISON  IS 
BETWEEN  A  SUM-OF-PAIRS  REPRESENTATION  AND  THESE 
THREE-BODY  FITS  TO  SECOND  ORDER 


In  Figure  13  the  influence  of  four-body  terms  (F^  )  is  shown. 

In  this  fit,  the  three-body  coefficients  of  fit  are  used,  and  only 
the  four-body  coefficients  are  adjusted  by  regression.  It  can  be  seen 
that  the  four-body  contribution  is  negligible  except  between  0°  and  30°, 
and  there  it  is  small.  This  is  the  case  for  all  other  cuts  that  were 
examined,  and  is  the  result  of  having  obtained  a  good  three-body  fit. 

In  Figure  14,  the  incoming  0-atom  is  fixed  at  various  distances 
from  the  molecule  cen ter-of-mass ,  and  the  molecule  is  bent  to  a  linear 
configuration.  The  four-body  fit  is  shown,  and  a  small  barrier  is  seen 
to  occur  for  the  3.0  aQ  curve.  This  is  an  artifact  produced  by  not  having 
a  sufficient  number  of  ab  initio  points.  In  Figure  15  the  differences  in 
the  fits  are  presented  for  a  symmetric  stretch  displacement.  The  three- 
body  curves  are  similar  in  shape,  but  the  four-body  curve  is  somewhat 
softer,  and  not  necessarily  more  accurate.  The  ab  initio  points  are  not 
shown  on  these  two  figures  because  the  geometries  used  do  not  exactly 
correspond  to  the  motion  described  in  these  figures.  Figure  16  presents 
a  yz  plane  view  of  the  potential  contours  for  the  0  +  H^O  surface, 
with  H2O  at  equilibrium.  The  surface  cut  is  seen  to  be  fairly  circular, 
only  skewed  somewhat  due  to  the  influence  of  the  hydrogen  atoms. 

Studies  of  many  cuts  through  the  various  surface  fits  demonstrate 
that  a  sufficiently  accurate  representation  of  the  Of  P)  +  HoO  PES 
is  obtained  (considering  the  limited  number  of  points  available)  using 
three-body  terms  through  second  order.  Quasiclassical  trajectories  are 
used  to  examine  differences  in  the  colli  si onal  excitation  cross-sections 
resulting  from  several  surfaces.  The  results  of  this  comparison  will  be 
given  in  the  last  section. 
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FIGURE  13.  SIMILAR  TO  FIGURE  11,  EXCEPT  FOR  THE  INCLUSION  OF  A  FOUR-BODY 
FIT  TO  SECOND  ORDER,  F,a .  THE  DISTANCE  TO  THE  CENTER-OF-MASS 
OF  THE  MOLECULE  IS  FIXED  HERE  AT  3.0  a  . 


FIGURE  14.  BENDING  FOR  THE  WATER  MOLECULAR  AS  A  FUNCTION 

OF  THE  APPROACH  OF  0(3P)  ALONG  THE  C2v  AXIS  TOWARD 
THE  OXYGEN  END  OF  THE  MOLECULE.  THE  CURVES  ARE 
ALL  FOR  THE  FOUR-BODY  FIT. 
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FIGURE  15.  A  COMPARISON  OF  THREE-  AND  FOUR-BODY  FITS  FOR  A 
SYMMETRIC  STRETCH,  WITH  0(3p)  3.0  a0  FROM  THE 
OXYGEN  AND  OF  THE  MOLECULE  ALONG  THE  C2V  AXIS. 


FIGURE  16.  POTENTIAL  CONTOURS  IN  THE  PLANE  OF  H20  FOR 
0  +  H2O.  THE  H2O  MOLECULE  IS  FIXED  AT 
EQUILIBRIUM.  CONTOUR  INCREMENT  IS  0.5  eV. 
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C.  Some  Aspects  of  the  Quasiclassical  Trajectory  Method  as 

Applied  to  Co*l  1  i si'onal  Excitation  of  H20 

Collisional  excitation  cross  sections  for  0  +  ^0  are 
determined  using  the  quasiclassical  trajectory  method.  Although 
the  basic  theory  of  trajectory  calculations  for  a  u>m- ;  6  tc/. 
collisions  was  developed  long  ago,  '  the  generalization  of  this 
to  collisions  involving  triatomics  and  higher  polyatomics  has  been 
investigated  in  detail  only  recently.  The  trajectory  integration 
itself  is  not  too  difficult  (although  it  can  involve  significant 
computer  time)  but  the  specification  of  the  initial  and  final  conditions 
of  the  polyatomic  collision  partners  requires  the  rather  difficult 

determination  of  their  semiclassical  vibration-rotation  eigenstates. 

( 27) 

If  the  molecular  internal  motion  is  quasiperiodic;  then  the  semi¬ 
classical  eigenstates  can  be  determined  by  finding  the  good  action 
variables  which  describe  vibration-rotation  motion,  and  equating  these 
actions  to  half  integral  multiples  of  h.  This  is  the  procedure  used 
for  each  trajectory  in  this  study,  both  for  the  initial  state  prepara¬ 
tion,  and  for  the  final  state  assignment. 

The  determination  of  the  good  actions  describing  vibration- 
rotation  motion  requires  the  solution  of  the  molecular  Hami lton-Jacobi 
equation  which  is  a  nonlinear  partial  differential  equation  in  3N-5  variables 
(including  rotation).  Even  for  N  =  3  (a  triatomic  molecule)  an  exact 
solution  to  this  equation  is  computationally  extremely  complex  and  is 
not  practical  for  collisional  applications.  Several  approximations  can 
be  used  to  simplify  this  treatment,  however,  including:  (1)  the  separation 
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of  vibration  from  rotation  (valid  in  the  limit  of  an  adequate  vibration- 

rotation  time  scale  separation);  and. (2)  the  use  of  classical  perturba- 

tion  theory  (in  2nd  and  3rd  order)  to  solve  the  three  dimensional 

Hamilton-Jacobi  equation  for  vibration  which  remains  after  the  separation 

from  rotation.  Details  of  both  the  separation  procedures  and  the 

(28  29) 

perturbation  theory  solution  are  discussed  elsewhere.  ’  For  the 

present  application,  the  validity  of  the  first  approximation  is  tested 

by  comparing  the  ensemble  average  change  in  translational  energy  with 

the  sum  of  the  vibrational  plus  rotational  energy  changes.  If  energy 

is  conserved  and  if  the  Hamiltonian  is  exactly  separable,  then  these  two 

energy  changes  should  have  the  same  magnitude  but  opposite  signs.  All  of 

our  calculations  indicate  that  this  holds  to  better  than  1%,  even  for 

highly  rotationally  excited  H2O.  The  accuracy  of  the  perturbation  theory 

method  is  tested  by  calculating  the  final  good  actions  from  a  given 

trajectory  at  several  different  end  points.  If  the  actions  are  truly 

constants  of  the  motion,  they  should  be  independent  of  end  point.  For 

H^O  with  zero  rotational  energy,  the  actions  oscillate  by  +  0.05R.  This 

is  much  larger  than  is  obta:ned  in  the  analogous  treatment  of  heavier 
(  ?r) 

molecules  like  C02>  but  is  consistent  with  expected  errors  due  to 
anharmonic  effects  in  H2O  and  is  probably  acceptable  for  calculating 
histogram  and  moment  probabilities.  When  the  H^O  molecule  has  rotational 
energy  equal  to  or  greater  than  the  vibrational  zero  point  energy,  the 
good  actions  oscillate  very  substantially  with  end  point  location  (by  as 
much  as  0.5h).  Obviously  such  large  oscillations  could  lead  to  inaccurate 
histogram  probabilities  and  moments.  Fortunately,  for  the  stretching 
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excitations,  the  amount  of  rotational  excitation  is  much  smaller  than 
the  zero  point  energy,  so  these  unphysical  oscillations  in  the  final 
actions  are  closer  to  0.05R  than  to  0.5R. 

One  very  effective  approach  for  partitioning  the  molecular 
Hamiltonian  to  simplify  vibration-rotation  interaction  effects 
involves  the  use  of  classical  rotation  sudden  approximations!^  These 
are  primarily  applicable  to  high-velocity  collisions,  and  have  been  quite 
effectively  applied  to  the  Li+  +  CO,,  system.  In  the  present  applica¬ 
tion,  both  sudden  and  fully  coupled  cross  sections  for  vibrational 
excitation  are  obtained. 

Usinc  either  the  sudden  or  fully  coupled  methods,  a  major 
problem  with  the  use  of  a  trajectory  method  to  calculate  state-to-state 
cross  sections  concerns  the  assignment  of  a  final  state  to  each  trajectory. 
One  well  studied  approach  is  the  classical  histogram  method.  In  this, 
the  final  good  actions  are  rounded  off  to  the  nearest  half  integer 
multiple  of  R,  and  the  trajectories  are  then  binned  by  the  final  states 
thus  defined.  This  is  known  to  work  well  when  many  states  are  populated 
in  the  excitation  process.  At  the  velocities  of  interest  in  this  study, 
however,  only  the  first  excited  symmetric  or  asymmetric  stretch  states 
are  found  to  be  excited.  This  makes  the  histogram  treatment  of  these 
excitations  somewhat  uncertain.  Moreover,  the  probability  of  stretch 
excitation  is  often  small  enough  in  this  system  (on  the  nrd er  of  0.001) 
that  the  accumulation  of  good  statistics  requires  the  integration  of  a 
very  large  number  of  trajectories  (many  thousands).  A  further  problem 
with  the  histogram  method  is  its  inability  to  determine  cross  sections 
for  transitions  which  are  classically  forbidden/  ; 


An  alternative  procedure  for  determining  cross  sections  (one 

which  circumvents  the  above  mentioned  problems)  involves  the  use  of 

moments  of  the  final  action  distribution  to  predict  the  final  state 

distribution.  In  essence,  this  "moment  method"  seeks  to  construct  the 

quantum  distribution  based  on  knowledge  of  a  few  classical  moments. 

Justification  for  this  procedure  is  based  on  the  approximate  equality 

of  low  order  classical  and  quantum  moments  for  some  systems,  and  on  the 

fact  that  it  is  possible  to  characterize  quantum  probability  distributions 

completely  in  terms  of  one  or  two  moments.  Although  moment  methods  have 

( 26  ^ 

often  been  studied  for  atom-diatom  systems;  '  with  comparisons  of  exact 
quantum  and  moment  method  probabilities  available,  no  analogous  study  of 
moment  methods  for  polyatomic  systems  have  been  done.  In  addition,  it  is 
not  particularly  clear  how  best  to  generalize  many  of  the  existing  moment 
methods  to  the  treatment  of  polyatomic  systems.  The  incorporation  of 
correlations  between  different  degrees  of  freedom  in  the  final  state 
distributions  has  not  been  considered,  for  example. 

In  view  of  the  desirability  of  a  moment  approach,  but  also 
the  uncertainty  in  the  application  of  existing  ones,  a  very  simple  method 
is  adopted  which  uses  moments  to  construct  and  extrapolate  the  classical 
action  distribution,  then  applies  the  histogram  method  to  the  moment 
generated  distribution  to  obtain  cross  sections.  Assuming  that  the 
exact  classical  action  distributions  can  be  accurately  represented  in 
terms  of  just  a  small  number  of  moments,  this  approach  will  generate 
moment  transition  probabilities  which  equal  the  histogram  ones  for 
classically  allowed  processes.  This  is  itself  quite  useful,  for  it 


provides  a  method  for  generating  cross  sections  for  low  probability 
allowed  processes  without  integrating  many  trajectories.  In  addition, 
it  provides  an  extrapolation  of  the  action  distribution  which  can  be 
used  to  estimate  weakly  forbidden  cross  sections. 

The  motivation  for  using  an  extrapolation  of  this  sort  is 
based  on  the  actual  appearance  of  the  action  distribution. 

When  summed  over  all  actions  except  the  one  of  interest,  this  distribution 
always  exhibits  an  exponential  decrease  with  action  on  both  sides  of  its 
peak  value.  For  any  finite  number  of  trajectories,  the  calculated 
distribution  will,  of  course,  truncate  when  the  distribution  decays 
to  a  low  enough  probability.  The  moment  method  smoothly  extrapolates 
th->  truncated  result  to  enable  calculations  of  low  probability  allowed 
and  forbidden  events.  For  allowed  processes,  the  method  has  obvious 
justification  while  for  forbidden  processes,  it  gives  a  smooth  extra¬ 
polation  of  the  allowed  results.  This  extrapolation  is  probably  inaccurate 
for  all  but  weakly  forbidden  processes,  but  since  the  transitions  of 
interest  in  this  work  are  always  either  allowed  or  weakly  forbidden  at  the 
velocities  of  interest,  the  extrapolation  is  small  and  presumably  accurate. 
For  most  of  the  moment  results  to  be  presented,  a  simple  exponential 
function  exp ( | N-N | )  is  found  to  describe  the  action  distribution 
adequately.  (Here  N  is  the  action,  N  the  average  of  N  and  y  is  a  decay 
parameter  which  can  be  related  to  the  second  moment  of  N).  Histogram 
probabilities  are  easily  determined  using  this  distribution  for  each 
mode,  and  the  total  (three  mode)  distribution  is  obtained  by  taking 
products  of  the  individual  mode  distributions. 
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In  this  study,  a  variable-step  fixed  order  Adams -Moulton 
method  is  used  to  integrate  the  equations  of  motion.  Maximum  separation 
is  fixed  at  10  aQ,  and  a  suitable  maximum  value  for  the  impact 
parameter  is  found  to  be  4.2  aQ.  The  initial  step-size  is  fixed  at 
five  atomic  seconds.  Sampling  of  impact  parameters  is  done  by 
randomizing  the  orbital  angular  momentum.  The  range  of  collision 
energies  studied  is  0.5  eV  to  2.5  eV,  and  usually  500  trajectories  are 
run  for  each  initial  translational  energy  and  rotational  state.  Initial 
rotational  energies  are  investigated  up  to  a  value  corresponding  to 
about  J=20. 
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D.  Cross  Sections  for  Collisional  Excitation 


Quasiclassical  trajectory  calculations  are  performed  on  a 
number  of  different  fits  to  the  ab  Initio  0 (  P )  +  HgO  ground  state 
potential  energy  surface,  using  full  three-dimensional  calculations  and 
a  classical  rotational  sudden  approximation.  The  use  of  several  fits 
allows  the  effect  of  variations  in  the  potential  surface  on  the 
computed  cross  sections  to  be  investigated.  In  all  cases  the  water 
molecule  is  in  its  ground  vibrational  state. 

In  Table  11  moment  cross  sections  are  compared  for  several 
different  fits,  all  involving  two-body  terms  to  second  order  This  particular 
translational  and  rotational  energy  are  chosen  because  results  are 
available  for  all  three  surfaces.  Considering  the  variations  in  the 
fits,  the  cross-sections  are  in  reasonable  agreement,  the  largest 


percent  error  being  for  the  asymmetric  stretch  between  fits  E.^  and 
F^i  but  still  less  than  a  factor  of  two.  Trajectories  are  not  run 
for  the  sum-of-pairs  representation  for  the  rotational  energy  in 


Table  11  ,  but  at  Er 


.16  the  001  cross  section  on  that  surface  is 


within  50%  of  the  cross  section  on  surface  E3a-  This  comparison  suggests 
that  errors  in  the  results  due  to  inaccuracies  in  the  fit  to  the  computed 


points  are  probably  less  than  a  factor  of  two. 

The  comparison  between  moment  cross  sections  and  histogram 
cross  sections  in  Table  12  shows  th-t  the  agreement  is  usually 
excellent  for  large  amounts  of  rotational  excitation.  For  collision 
energies  above  leV,  the  agreement  is  usually  better  than  a  factor  of  two. 
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TABLE 
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11.  COMPARISON  OF  MOMENT  CROSS-SECTIONS* 
FOR  SEVERAL  SURFACE  FITS 


Fit 

(100) 

(010) 

(001) 

E3a 

0.31 

0.99 

0.21 

F3b 

0.38 

0.73 

0.35 

F3c 

0.33 

0.66 

0.25 

ETRAN$  =  1,5  eV 
Erqt  =  0.08  eV 

Units  are  10"  ^  cmE 
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TABLE  12.  COMPARISON  OF  MOMENT  AND  HISTOGRAM  CROSS-SECTIONS* 


Er0T^v) 

Fn 

Sudden 

E3a 

0.08 

E,. 

3u 

0.16 

E3a 

0.40 

E3a 

0.77 

E2a 

010 

histogram 

0.65 

1.7 

1.4 

2.4 

3.4 

moment 

0.37 

0.99 

1.3 

2.3 

2.9 

001 

hi stogram 

0.06 

0.11 

0.22 

0.61 

0.90 

moment 

0.13 

0.21 

0.33 

0.58 

0.91 

TOO 

histogram 

0.22 

1.0 

1.1 

moment 

0.20 

0.31 

0.44 

0.87 

1.0 

E  =  1  5  eV 
tTRANS  ':)ev‘ 

- 1 6  2 

Cross-sections  in  10  cm  , 


I 

I 


This  is  not  the  case  below  1  eV,  at  least  in  part  due  to  the  poor 
statistics  for  the  histogram  results,  and  the  moment  results  are 
probably  more  reliable.  The  moment  results  are  generally  found  to  be 
reliable  only  for  one-quantum  transitions  once  the  histogram  probabilities 
become  large.  This  is  not  surprising  since  they  are  obtained  from  only 
the  first  two  moments. 

The  results  in  Table  12  also  demonstrate  the  almost  linear 
dependence  of  the  cross  sections  on  initial  rotational  energy.  These 
results  indicate  that  experimentally  determined  cross  sections  will  depend 
on  rotational  temperature,  and  that  experiments  performed  using  different 
(perhaps  unknown)  rotational  distributions  can  yield  considerably 
different  cross  sections.  The  sudden  results  are  found  in  general  to  be 
within  a  factor  of  two  of  tre  3-D  results  for  J  =  0. 

In  Figure  17  excitation  cross  sections  are  presented  for  an 
initial  rotational  energy  of  0.16  eV,  which  corresponds  to  about  J=  10. 

The  vertical  lines  through  some  of  the  points  give  lo  error  estimates  due 
to  Monte-Carlo  sampling.  Error  estimates  for  the  Oil  and  110  cross- 
sections  are  seen  to  overlap.  The  020  cross  section  is  not  shown,  for 
clarity,  since  it  almost  overlaps  the  001  cross  section  from  1-2  eV.  The 
200  and  002  cross  sections  are  estimated  from  the  moment  distributions, 
and  are  probably  low  by  a  factor  of  two.  Results  are  presented  at 
2.5  eV  only  for  the  100  and  010  cross  sections,  because  only  175  trajectories 
were  run  and  histogram  results  are  not  reliable  for  the  other  transitions. 
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Moment  cross  sections  (not  shown)  track  the  results  reasonably 
well.  For  example,  at  1.5  eV  the  moment  estimate  of  the  001  cross  section 
is  1.5  times  larger  than  the  histogram  result,  and  for  the  010  excitation 
the  histogram  result  is  1.2  times  larger  than  the  moment  result.  In  all 
cases  the  moment  results  are  smoother  as  a  function  of  translational 
energy. 

The  experimental  point  in  Figure  17  is  the  shock-tube  result  of 
Dunn  and  coworkers,  Reference  5,  for  the  excitation  of  the  001  asymmetric 
stretch.  Results  of  this  study  for  J  =  0  fall  almost  through  this  point, 
indicating  that  the  experimental  conditions  must  correspond  to  a  low 
rotational  temperature  to  be  consistent  with  these  calculations. 

Another  observation  from  the  present  results  is  pertinent  to 
the  interpretation  of  the  shock-tube  and  molecular-beam  experiments. 

Since  the  1 00  cross  section  is  greater  than  the  001  cross  section,  by 
about  a  factor  of  4,  and  since  the  emissions  from  the  two  states  overlap, 
even  though  the  100  state  has  a  much  longer  lifetime  (180  m  sec  compared 
to  17  m  sec),  the  larger  100  cross  section  might  contribute  about  40%  to 
the  observed  signal . 

In  summary,  the  results  obtained  in  this  study  indicate  that 
cross  sections  for  collisional  excitation  of  by  an  atomic  species  can 
be  obtained  to  within  about  a  factor  of  two  for  a  given  set  of  ab  initio 
points.  This  is  the  conclusion  reached  by  considering  different  fits  to 
the  points  as  well  as  two  different  techniques  for  generating  cross  sections. 


67 


In  summary,  the  results  obtained  in  this  study  indicate 
that  cross  sections  for  collislonal  excitation  of  H^O  by  an  atomic 
species  can  be  obtained  to  within  about  a  factor  of  two  for  a  yiven 
set  of  ab  initio  points.  This  is  the  conclusion  reached  by  considering 
different  fits  to  the  points  as  well  as  two  different  technique'  for 
generating  cross  sections. 

These  calculations  ignore  the  possibility  of  chemical  reaction. 

An  attempt  was  made  to  assess  the  importance  of  excited  electronic 

states  by  adjusting  the  sum  of  pairs  fit  to  correspond  with  the  cwo 

excited-state  surfaces  where  points  are  available.  Asymmetric  stretch 

2 

cross  sections  for  the  B2  state  (C2y)  are  similar  to  those  of  the 

2 

ground  stcte,  and  cross  sections  for  a  model  A^  surface  (C2y)  are 
about  30%  larger  than  the  ground  state  at  2eV.  Certainly  further 
work  on  this  system  must  involve  a  search  for  reactive  pathways  in  the 
potential  energy  surface.  Applications  of  the  surface  fitting  procedure 
and  the  dynamical  calculations  to  a  reactive  surface  would  be  tedious 
but  tractable. 


IV.  0(3P)  +  CCU1!*)  COLLISION 
A.  Potential  Energy  Surface 

The  collision  between  0(  P)  and  C02(  r  )  differs  in  several 
important  respects  from  that  obser  ved  with  H^O.  Primary  among  these 
are  that  the  CO  bond  is  much  stronger  than  an  OH  bond.  For  example, 

o 

the  stretching  force  constant  for  CO2  is  16  md/A  compared  to  8.4  in 
H^O.  From  another  viewpoint,  it  requires  more  than  twice  as  much 
energy  to  break  the  CO  bond  as  an  OH  bond.  These  facts  suggest  that 
COg  should  be  harder  for  an  oxygen  atom  to  vibrational ly  excite  than 
is  H^O.  To  the  contrary,  the  effective  radius  of  C02  is  substantially 
larger  than  in  H20.  However,  this  factor  contributes  only  to  the  elastic 
cross  sections.  Since  C02  is  a  linear  molecule,  this  surface  benefits 
from  higher  symmetry  than  the  H20  case,  enabling  us  to  obtain  more 
information  from  each  computed  point.  Except  for  approaches  along  the 
axis,  the  0(  P)  +  C02(  e  )  correlation  diagram  is  similar  to  that 
of  H20.  This  is  illustrated  in  Figure  18.  Figure  19  shows  the  orienta¬ 
tion  of  C02  in  the  yz  plane  and  the  distribution  of  ab  initio  points 
on  the  surface. 

When  the  oxygen  atom  approaches  along  the  z-axis  (i.e.,  the 

axis),  the  possible  electronic  states  are  a  1  and  a  doubly  degenerate 
3 

n.  Just  as  in  H^O,  an  approach  along  the  Cr,  axis,  which  is  the  y-axis 

3  3  3 

in  Figure  19,  generates  the  three  states,  ,  B2>  and  A2 .  However, 

unlike  H20,  when  C02  is  stretched,  symmetrically  or  asymmetrically  (but 
not  bent),  the  values  computed  at  the  points  indicated  on  the  y-axis  define  a 
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0  +  COg 

ob  initio  Point  Distribution 
Rc,0  =  2.29I5B 
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FIGURE  19.  DISTRIBUTION  OF  AB  INITIO  POINTS  IN  THE  PLANE 
OF  THE  C02  MOLECULE,  FOR  EQUILIBRIUM  C02 
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locus  of  circles  about  the  carbon  atom.  When  a  single  plane  of 
symmetry  is  retained,  such  as  by  any  approach  in  one  of  the  planes 
of  the  CO2  molecule,  the  controlling  symmetry  is  Cs,  providing  a 
ground  A'  and  two  excited  A"  states. 

As  justified  in  our  study  of  the  0(  P)  +  HgCK  A^ )  surface, 
there  seems  to  be  little  reason  to  determine  the  0(  P)  +  C02 (  I  ) 
surface  using  the  much  more  expensive  DZP  basis  set.  Consequent iy, 
we  exclusively  employ  a  double-zeta  ( DZ }  basis  in  this  part  of  the 
project.  Such  a  basis  is  largely  incapable  of  showing  much  of  the 
very  slight  Van  der  Waals  region,  which  should  occur  about  7-8  bohrs 
in  this  system.  In  the  DZP  studies  of  0(  P)  +  H^O,  this  region  is 
somewhat  described,  although  even  a  DZP  basis  is  too  small  to  provide 
an  accurate  description.  However,  a  fit  could  introduce  the  slight 
Van  der  Waals  well  for  0  +  CO^  from  appealing  to  experiment.  In  the 
current  studies,  this  region  is  found  to  have  such  a  small  effect  on 
the  cross  sections  that  it  is  not  considered  an  important  element  in  the 
surface  in  any  event.  Other  than  the  Van  der  Waals  region,  the  SDQ-MBPT(4) 
DZ  calculations  are  found  to  parallel  the  DZP  calculations  on  0  +  H2O. 

The  predicted  equilibrium  structure  for  C02  gives  a  bond 
length  of  2.2915  bohrs  which  differs  oy  about  0.08  bohrs  from  the 
experimental  value.  The  separated  limit  is  -262.74654  hartrees. 

With  reference  to  the  distribution  of  ab  initio  points  for 
CO2  at  equilibrium,  shown  in  Figure  19,  the  computed  values  are 
presented  in  Table  16.  Figure  20  illustrates  the  four  curves  for 
different  angles  of  approach.  Like  the  H20  +  0(  P)  surface,  the 
repulsive  nature  of  the  curves  is  evident.  In  H2O,  the  surface 
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is  almost  spherical,  while  for  CO^  Figure  10  emphasizes  the 
cylindrical  character  that  one  would  expect. 

The  steepest  repulsive  shoulder  corresponds  to  a 
collision  along  the  0^  axis,  rising  0.19  a.u.  (i.e.,  5 eV)  from 
7  to  5  bohrs.  This  may  be  compared  with  the  y-axis  approach  toward 
H20,  illustrated  in  the  bottom  curve  of  Figure  5,  which  is  seen  to 
be  less  abrupt  from  6  to  4  bohrs.  This  shows  that  the  oxygen  atom 
may  approach  the  smaller  molecule  more  closely,  but  then  increases 
abruptly  by  0.35  a.u.  in  the  one  bohr  from  3.5  to  2.5  bohrs.  Since  the 
H-atoms  are  displaced  from  the  y-axis  in  H^O  by  about  38°,  another  pertinent 
comparison  is  between  the  ft  =  30°  curve  and  Figure  5.  In  this  curve, 
the  approach  is  less  abrupt,  changing  0.24  a.u.  from  7  to  4  bohrs,  compared 
to  almost  twice  that  in  H20,  between  5.5  and  2.5  bohrs.  As  ft  increases 
to  90°,  the  curves  are  shaped  more  like  those  in  the  H.,0  collision,  but 
still  remain  a  little  less  steep,  beginning  their  ascent  in  the  vicinity 
of  5  bohrs  and  rising  to  0.2  a.u.  in  2.5  bohrs.  The  slopes  of  the 
four  curves  of  Figure  20  in  these  steep  regions  are  all  quite  similar. 

Through  the  ability  of  the  oxygen  atom  to  approach  H20 
more  closely  than  C02>  a  more  drastic  increase  in  the  slope  occurs, 
making  the  H20  collision  harder  at  short  distances,  while  softer  in 
the  intermediate  range.  The  latter  is  due  in  some  degree  to  the 
influence  of  the  dipole  moment  of  H^O. 

In  Tables  14  to  16,  the  computed  points  for  symmetric  and 
asymmetric  stretch  displacements  and  a  few  points  for  bent  C02  are 
reported.  For  a  symmetric  displacement  the  values  on  the  z-axis  are 
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TABLE  14.  0(3p)  +  C02( 1 )  GROUND  STATE  SURFACE  FOR  C02  IN  A  SYMMETRIC  STRETCH 

(Rrn  =  2.35  b,  9nrn  =  180°) 


9  m  ■  » w 


TABLE  15.  0(JP)  +  C02  (  I*)  GROUND  STATE  SURFACE  FOR 

C02  IN  AN  ASYMMETRIC  STRETCH  MODE. 

[O^O  ,0  ,  -2.33242b)  C(0.,  0.,  +  0.10908b), 
09(0.,  0.,  +  2.25060b)] 


Points 


.0 

90 

- 

262.26705 

-262.66806 

80 

90 

- 

262.26211 

-262.66304 

90 

90 

- 

262.26947 

-262.66471 

90 

90 

- 

262.30853 

45 

90 

- 

262.26587 

-262.66962 

45 

90 

- 

262.30759 

-262.70656 

35 

90 

- 

262.26729 

-262.67232 

35 

90 

- 

262.30822 

Angles  0  =  45 


TABLE  16.  0(3P)  +  C02(1lg)  GROUND  STATE  SURFACE  FOR 

C02  IN  A  BENT  CONFIGURATION  [0]  (0-j ,  0.238628b, 
-2.255210b),  C(0.,  -0.159031b,  0.), 

02(0.,  0.238628b,  +  2.25521b)] 


Poi nts 


r 

9 

* 

SCF 

SDQ-MBPT(4) 

2.7 

90 

90 

-262.06239 

-262.47364 

3.0 

90 

90 

-262.15400 

-262.56164 

5.3 

0 

90 

-262.23997 

-262.64853 

5.5 

0 

90 

-262.27269 

0262.67856 
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slightly  higher  partially  reflecting  a  slight  increase  in  energy  due  to 
CO2  being  displaced  from  its  computed  equilibrium  value,  and  partially 
due  to  the  oxygen  atom  and  molecule  being  somewhat  closer.  The  same  trend 
applies  in  the  case  of  the  30°  and  60°  approaches,  although  at  n  =  90°, 
the  symmetrically  displaced  terminal  atoms  serve  to  slightly  facilitate 
the  nearer  approach  of  the  oxygen  atom. 

In  the  asymmetric  case,  results  are  obtained  that  largely 
depend  upon  which  side  of  the  molecule  is  nearest  to  the  colliding 
oxygen  atom.  A  difference  of  .005  hartrees  is  observed  on  the  2-axis 
at  5.5  bohrs,  and  somewhat  less  for  the  approaches  at  45°  and  135°,  slightly 
favoring  the  latter.  In  the  asymetric  stretch  configuration,  the  energies 
on  the  z-axis  are  both  higher  than  their  values  with  C02  at  equilibrium. 

Just  as  in  0(3P),  the  ground  state  surface  for  0(3P)  +  CO^r*), 
shows  a  comparatively  simple  repulsive  nature,  which  can  be  conveniently 
fit  to  an  analytic  fonn  as  discussed  in  Section  IV.  B.  The  excited  states 
that  also  correlate  with  an  Q(  P)  atom  at  the  separated  limit  are  not 
considered.  In  the  case  of  H^O,  these  are  found  to  make  very  little 
difference  in  the  cross  sections.  The  program  plan  specifies  studying 
collisions  only  on  the  ground  state  surface;  and,  because  the  cross  sections 
obtained  in  Section  IV. D.  are  so  small,  a  more  detailed  study  incorporating 
the  excited  spates  would  be  very  unlikely  to  gain  more  than  a  factor  of  3 
or  so,  while  at  least  an  increase  of  one  or  two  orders  of  magnitude  would 
be  necessary  to  make  the  direct  excitation  of  the  CC^lOOl)  mode  (as  opposed 
to  indirect  excitation)  an  important  contributor  to  the  infrared  plume 
signature. 


79 


B.  The  Analytical  Representation  of  the 
0(3P)  +  COp  Potential  Energy  Surface 

The  PES  for  0(  P)  +  CO^  is  represented  in  the  same  general 
form  as  the  0(  P)  +  H20  surface,  namely  as  the  sum  of  intramolecular  (C02) 
and  intermolecular  terms.  The  intramolecular  potential  for  CO2  is  represented 
as  a  many-body  expansion  similar  to  H20.  There  are  two  significant  differences, 
however,  the  first  being  the  use  of  an  experimental  CO^  force  field, 
rather  than  a  theoretical  force  field.  The  second  is  the  use  of  internal 
symmetry  coordinates  instead  of  bond  displacements  for  the  CO^  polynomial. 

(00) 

The  experimental  force  field  of  Kuchitsu  and  Morino  is  used  since, 
unlike  H^O,  no  theoretical  determination  of  the  C02  force  field  of  sufficient 
accuracy  exists. 

The  form  of  the  C02  potential  is 

VC0^R1’  R2’  V  =  WV  +  VC0^R3^  +  V02^V  +  W  R2’  V 

where  R-j  and  are  the  two  C-0  bond  distances  and  R2  is  the  0-0  bond 

1  3  - 

distance.  VCQ  and  VQ  are  the  C0(  l)  and  02(  1  )  diatomic  potentials. 

g 

is  a  three-body  term  and  has  the  form 

V3  =  A(1  -  tanh  -,-jS^/2)(l  -  tanh  y2s|/2 ) ( 1  -  tanh  y2S2/2)P(S-|,  S2,  S3) 

where  A,  ,  yn,  and  73  are  constants,  and  P  is  a  polynomial  in  the 
internal  coordinates.  These  are  defined  as 
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(R*  -  R?)/(2R2) 


(1  r2  +  1r2  .  1  R2  .  s2}1/2 

l2  R1  2  R3  T  2  V 


s3  =  r2  ‘  *2 


where  is  the  molecular  equilibrium  value  of  R2>  Note  that  vco2  should 

be  symmetric  with  respect  to  S -j S 3  and  S2+-S2.  For  this  reason  the 

2  2 

switching  functions  in  V3  are  chosen  to  be  functions  of  S-j  and  S,,  rather 

than  S-j  and  S2-  The  use  of  S-| ,  »  and  rather  than  R^,  and  R^»  as 

is  done  for  H20,  is  necessary  because  the  polynomial  expansion  of  V3  in 

the  latter  variables  provides  an  ill-behaved  representation  of  the  quartic 

force  field  for  linear  molecules. 

The  parameters  for  V'cq  are  given  in  Table  17.  In  Tablel8 

the  quartic  force  constants  obtained  from  Vr«  are  compared  with  those  of 

1/U2 

Kuchitsu  and  Morino,  and  the  agreement  is  seen  to  be  excellent.  The 
present  fit  should  provide  a  realistic  extrapolation  from  the  accurate 
representations  around  equilibrium  to  the  separated  diatomic  fragments. 

Figure  21  shows  potential  contours  for  stretching  displacements  of  R-|  and  R3. 

.  3 

The  0  +  C02  mtermolecular  potential  for  the  lowest  A"  surface 

is  represented  as  a  sum-of-pairs  potential,  similar  to  the  E2g  fit  for 

H20.  Time  did  not  permit  including  higher  order  terms,  but  experience  with 

this  level  of  fitting  indicates  that  a  reasonable  representation  of  a 

four-atom  nonreactive  surface  is  obtained  in  this  way.  Excitation  cross 
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TABLE  17.  PARAMETERS  IN  V, 


(a) 


co2 


A.  Two  Body  Terms 

Vxy  =  'DXV(1  +  alxyARxy  +  a2xy  ARxy  +  a3xy  ARxy)e  ***** 


2L 

5<y 

alxy 

a2xy 

a3xy 

R° 

CO 

11.2245 

3.7102 

1 .5944 

1 .5422 

1.1282 

00 

5.2130 

5.4454 

7.7859 

5.9951 

1.2078 

B.  Three  Body  Terms 


A  =  3.9909 
Y-,  =  3.571068 

y2  =  3.571068 
y3  =  5.4454 
=  2.32303 
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TABLE  17.  (Continued) 


Polynomial  Term  Coefficient 


1 

1 

S3 

1  .5678 

bl 

1  .5871 

s2 

b2 

2.8483 

S3 

5.6727 

slS 

2.8466 

S22S3 

3.1655 

S3 

7.5278 

-4.2510 

s  2S2 

bl  b2 

3.9974 

s  2S2 
bl  ^3 

6.8399 

s4 

^2 

3.1893 

S  2$2 
^3 

18.391 

S4 

8.8036 

O 


Units  are  such  that  unprinted  distances  in  A  give  energies  in  eV. 


TABLE  18.  QUARTIC  FORCE  FIELD  PARAMETERS 
FOR  C02  POTENTIAL  FIT ( 3 ) 


Present  Fit 

Kuchitsu  and  Morino 

W1 

1354.8 

1354.8 

“2 

672.7 

673.2 

w3 

2395.7 

2395.7 

kni 

45.0 

45.0 

k122 

-73.9 

-73.8 

kl  33 

252.5 

252.5 

knn 

1.2 

1.2 

k3333 

6.6 

6.6 

kl  1 22 

-9.1 

-9.1 

kl  1 33 

20.3 

20.3 

k2233 

-27.7 

-27.8 

^  All  parameters  are  in  cm"\  For  definitions 
of  terms  used,  see  Reference  25. 
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0-0 


sections  for  0  +  H20  on  a  sum-of-pairs  surface  differ  by  no  more  tnan 
a  factor  of  two  from  cross  sections  on  a  more  accurate  fit  containing 
three- body  terms. 

Denoting  the  0  atom  to  C  distance  as  R^,  and  the  two  0  atom 
to  0  distances  as  R5  and  Rg,  the  expression  for  the  intermolecular 
potential  is 

-aR,  -BRc  -SRg 

voco2  =  Ae  +  Be  +  Ce 

The  parameters  in  this  expression  are  A  =  286.97  eV,  B  =  2195.9  eV, 
a  =  1.5563  a'1 ,  and  B  =  2.2250  a'1. 


C.  Quasiclassical  Trajectory  Calculations  as  Applied 
to  Collisional  Excitation  of  COn 

Both  full  3-D  and  rotational  sudden  calculations  are  made 
for  0(JP)  +  COg  as  f°r  0(3P)  +  H2O.  The  procedure  for  integrating 
trajectories  is  similar,  but  there  is  one  difference  in  treating  the 
initial  and  final  conditions.  For  CC^,  an  approximation  is  introduced 
that  involves  the  omission  of  a  certain  Fermi  resonant  term,  one  that 
couples  the  bend  and  symmetric  stretch  modes  in  the  Fourier  representation 
of  the  good  action  variables.  This  term  cannot  be  treated  using  the 
presently  available  nondegenerate  perturbation  method  for  defining  action- 
angle  variables.  Severa1  numerical  tests  indicate  that  this  term  has  no 
influence  on  excitation  of  the  asymmetric  stretch  mode.  Note  that  this 
approximation  is  only  introduced  to  allow  the  use  of  perturbation  theory 
to  determine  the  good  actions  used  in  defining  quantum  states,  and  that 
the  full  interaction  potential  is  used  to  integrate  the  trajectories. 

As  discussed  in  the  application  to  HjO,  a  problem  arises  with  the 
use  of  the  histogram  method  in  defining  state-to-state  cross  sections  when 
very  few  or  no  trajectories  finish  in  the  desired  bin.  This  is  often  true 
for  the  C02  asymmetric  stretch  excitation,  since  rarely  more  than  a  few 
trajectories  (out  of  1000  to  2000) lead  to  asymmetric  stretch  excited 
states.  To  extrapolate  and  interpolate  cross  sections  at  velocities 
where  poor  statistics  make  the  historgram  results  unreliable,  a  moment 
analysis  can  be  employed,  and  checked  against  histogram  results  where 
they  are  available. 


°7 

u/ 


In  0  +  H  0,  an  exponential  expression  is  found  to  conveniently 
2 

represent  the  final  state  distributions.  In  the  case  of  0  +  CC^.  this 
is  not  as  satisfactory,  but  an  expression 

?  x 

0  a  <aN^> 

seems  to  be  suitable.  The  value  X  =  1.3  accurately  represents  the 
nonzero  histogram  results  and  is  consequently  used  to  determine  the 
interpolated  cross  sections. 


D,  Cross  Sections  for  Collisional  Excitation 
of  CQo  by  o(3p)  Atom  Impact 

Figure  22  and  Table  19  presents  the  theoretical  cross  sections  For 
excitation  of  the  asymmetric  stretch  mode  of  CO2 .  The  result  levelled 
Q ( 001 )  i s  for  the  direct  process  000->001 ,  while  the  label  Q(Nfll)  results  from 
summation  over  the  symmetric  stretch  quantum  number  N  and  the  bend  quantum 
Number  N  for  those  transitions  leading  to  one  quantum  of  excitation  in 
the  asymmetric  stretch  mode.  Since  infrared  fluorescence  between  (NN'l) 
combination  states  and  (NN  0}  significantly  overlap  the  001  ->000  emission, 
it  is  likely  that  many  experiments  are  actually  measuring  the  total  cross 
section  Q ( NN1  )  rather  than  Q ( 001 ) .  At  velocities  of  6  km/sec  or  greater, 
this  study  shows  that  combination  state  excitation  is  substantially  more 
proDable  than  direct  excitation  of  the  001  state,  so  that  p(NN  1  )  is  much 
larger  than  Q ( 00 1 ) . 

The  results  in  Table  19  suggest  that  Q ( 001 )  is  immeasurably 
small  except  at  very  high  velocities.  At  6  km/sec,  Q ( 001 )  has  a  value 
from  this  study  of  3  x  10  cm  ,  while  Q(NN'l)  has  a  much  larger  value 

TO  0 

of  1  x  10  cm  .  Statistical  uncertainty  in  these  results  is  +  50  peuent. 

The  effects  of  CO2  rotational  excitation  on  both  Q(001)  and 
Q(NN'l)  are  summarized  in  Table  20.  A  slow  increase  in  these  cross 
sections  with  increasing  J  is  found  up  to  high  rotational  energy.  Although 
this  result  is  different  than  found  for  0  +  H^O ,  it  is  consistent  with 
tne  greater  time  scale  separation  between  rotation  and  vibration  in  C0o 
relative  to  H^O  due  to  the  large  CO2  moment  of  inertia. 
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TABLE  20.  EFFECT  OF  C02  DOTATION  ON  VIBRATIONAL 
EXCITATION  CF.OSS  SECTIONS(a) 


ER0T(eV) 

0 

.025 

j(co2) 

0 

23 

42 

No.  trajectories 

2000 

1000 

1000 

MAX  l 

550 

550 

550 

Q(001) 

2.5  x  10‘19 

.19(b) 

2.8  x  10  Iy 

_1< 

3.9  x  ’0 

Q(NN'l) 

8.8  x  10‘18 

..17(b) 

1.0  x  10  " 

1.4  x  10"1, 

a  Cross  sections  are 

2 

in  cm  .  Relative  velocity  is  8.1 

km/sec. 

b  Result  of  scaling 

calculation. 

Where  histogram  values  were  also 

obtained,  the  scaling  result  agrees  with  the  histogram  result 
within  the  statistical  uncertainty  of  the  latter. 


V.  DISCUSSION  OF  RESULTS  AND  SUGGESTIONS  FOR  FURTHER  WORK 


This  report  presents  results  of  a  first-principle  theoretical 
determination  of  cross  sections  for  vibrational  excitation  of  and  CO^ 
from  collisions  with  0(  P)  atoms.  Correlated  potential  energy  surfaces 
are  obtained  with  many-body  perturbation  theory  and  fit  to  analytical 
expressions.  Tne  analytical  fits  are  used  in  quasiclassical  trajectory 
calculations  to  determine  vibrational  excitation  cross  sections  averaged 
over  final  rotational  states.  Classical  perturbation  theory  is  used  to 
solve  the  Hami 1  ton- Jacobi  equations  to  obtain  good  actions  for  the  semi- 
classical  final  state  analysis. 

A  goal  in  this  work  is  the  prediction  of  cross  sections  to 
within  a  "factor  of  two"  of  their  accurate  values.  To  this  end  the  effects 
of  errors  in  tne  results  due  to  approximations  in  each  phase  of  tne 
calculation  are  investigated  by  running  trajectories  on  various  fits  to 
the  potential  surfaces.  This  proced"re  uses  fits  that  differ  from  one 
another  by  more  then  the  errors  In  an  ab  initio  surface,  and  also  fits 
that  reflect  the  differences  between  the  ground  state  surface  and  the 
first  two  excited  states  for  ^0.  In  most  cases  the  cross  sections 
are  changed  by  less  than  a  factor  of  two.  One  exception  occurs  when 
a  quartic  force  field  is  used  for  the  H^O  potential,  in  wnich  case  the 
cross  sections  differ  from  a  more  realistic  fit  by  an  order  of  magnitude. 
Convergence  to  the  correct  potential  is  assured  by  comparisons  of 
different  levels  of  correlation  and  by  using  different  basis  sets. 


Theoretical  studies  such  as  this  are  not  only  of  value  because 


of  the  large  number  of  excitations  for  which  results  are  obtained,  but 
as  an  aid  to  the  experimentalist  as  well.  In  addition  to  providing  results 
for  transitions  that  cannot  be  determined  by  experiment,  results  often 
can  help  in  the  interpretation  of  the  experiments.  The  dependence  of  the 
cross  sections  on  rotational  energy  discussed  in  the  text  indicates  that 
reagent  rotational  temperature  will  affect  the  vibrational  energy  transfer. 

This  might  be  important  in  interpreting  shock-tube  results.  Also  a  larger 

cross  section  for  100  than  001  in  H?0  can  affect  the  experimental  interpretation. 

The  results  for  CCL  show  that  excitation  of  Hfil  states  is  more 
probable  by  up  to  two  orders  of  magnitude  than  direct  excitation  of  the 
001  state.  As  discussed  in  the  text,  radiation  from  NNl  states  is  likely 
to  obscure  the  001  radiation,  and  this  can  affect  cross  sections  obtained 
in  some  experiments.  If  radiation  from  001  is  indeed  observed,  then  a 
multiple  collision  environment  probably  exists  and  a  study  of  collisional 
deactivation  of  excited  vibrational  states  to  001  is  necessary.  It  should 
be  noted  that  final  state  rotational  distributions  and  angular  distributions 
can  be  obt: ''  led  from  the  theoretical  approach,  and  would  likely  lead  to 
an  increased  understanding  of  these  collision  processes. 

Just  as  there  are  many  difficulties  in  complex  beam  and  shock 
tube  experiments,  there  are  characteristic  problems  in  the  theoretical 
calculations  as  well.  These  problems  are  discussed  more  fully  throughout 
this  report,  but  it  is  useful  to  summarize  the  possible  errors  here  in 
order  to  assess  the  reliability  of  the  computed  results  and  to  suggest 
possible  future  refinements  of  the  calculations. 
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Initially,  it  should  oe  emphasized  that  all  results  reported 
here  are  first  principle,  except  that  empirical  information  is  used 
for  the  CO2  force  field.  The  information  with  which  we  start  is  simply 
the  nuclear  charges  and  the  numbers  of  electrons  of  the  atoms  involved. 

The  electrostatic  hamiltonian  constructed  from  this  information  includes 
the  kinetic  energy  of  the  electrons,  the  electron-nuclear  attraction, 
and  the  interelectron  repulsions.  Armed  with  such  a  hamiltonian 
we  attempt  to  solve  Schrodinger's  equation  to  define  the  PES. 

The  method  of  solution  is  approximate  in  the  sense  that  a  finite 
number  of  basis  functions  are  used  and  all  possible  interelectronic  repulsive 
effects  (called  correlation)  are  not  included.  However,  demonstrated 
experience  tells  us  that  the  level  at  which  we  are  solving  Schrodinger’s 
eauation  is  sufficient  for  the  problem  at  hand.  We  also  find  only 
insignificant  differences  in  the  PES  when  chanqinq  from  a  small  number  of 
basis  functions  (i.e.,  the  double-zeta  level,  DZ)  to  a  larger  number 
(i.e.,  double-zeta  plust  polarization,  DZP).  Similarly,  only  small  changes 
occur  between  the  uncorrelated  self-consistent  field  results  and  our 
SD0-MBPT(4)  model,  which  essentially  adds  all  the  correlation  effects 
due  to  single-,  double-,  and  quadruple-excitations  into  the  problem. 

This  tells  us  that  an  improvement  in  basis  set  or  in  a  hiqher  inclusion 
of  correlation  effects  is  not  going  to  make  any  important  difference  in 
the  PES,  so  Schrodinger's  equation  is  essentially  solved  for  the  ground 
state  PES. 


95 


The  next  stage  of  the  calculation  involves  the  analytic  fit 
to  the  computed  ab  initio  points  whose  accuracy  seems  sufficient.  The 
accuracy  of  the  cross  sections  is  therefore  primarily  a  function  of  the 
quality  of  the  fit  to  the  computed  points,  if  the  fit  is  good  most  places, 
but  introduces  wells  or  barriers  that  should  not  be  present,  cross 
sections  could  be  inaccurate.  This  is  a  difficult  part  of  the  problem 
because  as  yet,  there  is  no  general  procedure  for  fitting  a  four-atom 
surface. 

As  a  result  there  are  many  approaches,  and  building  upon  one 
of  these  that  has  been  found  to  be  successful,  a  satisfactory  procedure 
has  been  developed  during  the  course  of  this  work.  We  have  been  greatly  helDed 
in  this  by  the  simple  repulsive  nature  of  the  computed  surfaces.  About 
tne  only  test  of  the  fit  that  is  possible,  other  than  simply  making  several 
plots  that  compare  the  fit  to  the  computed  points,  is  to  attempt  some 
"sensitivity"  analysis.  Essentially,  this  means  compute  the  cross-sections 
from  trajectory  calculations  for  a  given  fit  to  the  surface,  then  modify 
several  of  the  parameters  in  the  fit  and  repeat  the  cross  section  calculations. 
If  tnis  is  done  *or  a  fairly  wide  range  of  possible  fitting  parameters, 
it  is  possible  to  conclude  that  the  cross-sections  change  within  a  specified 
tolerance.  This  is  done  for  the  H^O  cross  sections,  concluding  that  the 
cross  sections  are  insensitive  to  within  a  factor  of  2  for  substantial 
changes  in  the  fit  parameters.  This  provides  a  measurable  degree  of 
confidence  in  the  predicted  results.  Due  to  the  use  of  a  simpler  fit  for 
the  0  +  C09  surface,  the  cross  sections  are  not  as  reliable  as  for  0  +  H^O. 
Future  work  should  involve  obtaining  a  better  fit. 


The  last  element  in  the  computation,  the  actual  determination 
of  the  cross  sections,  also  involves  approximations.  Recognizing  that 
the  collision  velocities  of  interest  in  this  work  are  relatively  high, 
it  is  possible  to  use  quasiclassical  trajectory  methods  to  obtain  the 
cross  sections.  This  means  using  classical  mechanics  to  describe  the 
actual  colliding  particles,  but  combining  these  techniques  with  good 
action-angle  variables  that  permit  the  trajectories  to  be  reliably 
resolved  into  specified  vibrational-rotational  states.  When  the  number 
of  trajectories  binned  into  a  vibrational-rotational  state  is  statisti¬ 
cally  inadequate,  tien  either  an  exorbitant  number  of  trajectories  must 
be  computed  to  obtain  statistically  meaningful  results,  at  great  expense, 
or  the  histogram  method  needs  to  be  augmented  with  moment  analysis,  as 
described  in  this  report.  The  application  of  moment  methods  to  polyatomics 
is  new  with  this  work,  and  its  reliability  for  different  types  of  excitation 
nas  not  been  investigated.  Care  is  taken  in  this  study  to  obtain  pro¬ 
cedures  that  are  reliable  for  the  asynmetric  stretch  excitations. 

Since  there  are  many  possible  rotational  states  for  each 
v’bration,  several  thousand  trajectories  would  be  required  for  statistically 
meaningful  results  for  rotational  distributions.  The  objective  of  this 
work,  instead,  is  to  provide  vibrational  cross  sections  averaged  over 
the  rotational  states.  For  the  case  of  a  very  small  collision-induced 
vibrational  excitation,  such  as  the  (001)  cross  section  in  CO,,  statistical 
uncertainties  can  be  high,  but  recognizing  the  small  magnitude  of  these 
quantities  is  sufficient  to  discount  the  importance  of  direct  excitation 
of  this  state  in  the  p.ume  signature. 
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The  other  approximation  in  the  actual  cross  section  determi¬ 


nation  involves  the  perturbation  approach,  required  for  efficiency,  that 
is  used  to  resolve  the  trajectories  into  the  various  quantum  states.  For 
very  high  internal  excitation,  the  perturbation  theory  does  not  converge. 
Also,  when  the  amount  of  rotational  energy  is  equivalent  to  vibrational 
excitation,  coriolis  coupling  tends  to  produce  large  temporal  oscillations 
in  the  actions.  These  problems  do  not  affect  the  present  results. 

This  discussion  describes  some  o*  the  significant  approximations 
made  in  this  work,  all  of  which  we  believe  are  sufficiently  accurate  for 
tnis  study.  However,  there  are  other  apDroximations  which  are  more 
difficult  to  discount  without  additional  work.  Since  all  of  these  are 
outside  the  scope  of  the  current  project,  and  since  they  mostly  pertain 
to  nigner-order  effects  on  the  cross  sections,  they  have  only  been 
considered  superficially.  A  more  thorough  investigation  is  warranted. 

Since  an  0(~P!  atom  colliding  with  a  molecule  in  a  totally 
symmetric  ground  state  produces  three  electronic  PES,  the  cross  sections 
for  the  collision  are  a  function  of  all  these  states.  One  computes  a 
cross  section  for  each  possible  electronic  PES,  aid  the  total,  observed 
state-to-state  cross  sections  are  found  as  one-third  of  the  sum  of  the 
cross  sections.  The  difficulty  and  expense  of  determing  these  PES 
requires  that  this  first  effort  should  essentially  focus  only  on  the 
ground  state  as  it  does,  although  in  H?0,  the  excited  states  are  also 
considered  to  some  degree.  When  all  three  states  are  determined,  t.ie 
two  lowest  are  found  to  be  quite  similar,  with  the  third  somewhat  more 
repulsive.  By  adjusting  the  fit  of  the  surface  to  account  for  the 
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changes  in  shape  in  the  excited  surfaces,  an  estimate  of  the  change  in 
the  cress  section  due  to  excited  states  is  made,  and  found  to  he  on  the 
c*“der  of  30  percent  for  H^O. 

In  CO^ ,  only  the  ground  surface  is  considered.  If  there  are 
larger  changes  among  the  three  states  than  anticipated,  different,  probably 
somewhat  larger  cross  sections  would  be  found.  A  thorough  study  of  these 
excited  state  PES  would  require  an  approach  for  solving  Schrodinger ’ s 
equation  that  gives  excited  states  of  the  same  symmetry  as  the  lower 
states.  We  are  developing  many-body  approaches  for  this  problem,  Cl 
approaches  are  applicable,  and  some  MC-SCF  techniques,  but  this  is  still 
a  difficult  problem.  Indirect  pathways  for  producing  4.3  u  radiation, 
probably  involving  multiple  collisions,  should  also  be  investigated 

Another  area  for  a  more  thorough  investigation  pertains  to  the 

■3  i  ? 

possible  reactive  channel  in  0(VP)  +  0  (  )-*20H(T ) .  To  study  thir.,  wr  must 

focus  on  the  reactive  reqion  of  the  enerqy  surface  and  locate  a  transition 

2 

state.  It  appears  from  the  correlation  diagram,  that  the  0H(  \ )  molecule 
can  be  formed  on  both  the  ground  and  first  excited  PES,  again  placing 
empnasis  on  a  knowledge  of  the  excited  state.  The  sign  of  any  possible 
error  due  to  this  approximation  is  hard  to  assess.  The  effect  of  reaction 
on  the  nonreactive  flux  could  be  estimated  if  the  properties  of  the 
transition  state  were  known.  The  more  reaction  actually  occurs,  the  less 
flux  will  enter  the  vibrational  excitation  channels,  reducing  the  cross 
sections.  However,  if  the  surface  develops  a  pathway  leading  to  these 
products,  it  may  also  enhance  the  vibrational  cross  sections  by  providing 
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an  area  that  would  encourage  repeated  collisions.  Also,  the  possible 
infrared  emission  of  OH  itself  could  also  contribute  importantly  to 
the  plume  signature,  and  this  should  be  investigated  for  that  reason 
as  well . 

A  third  area  for  future  research  would  be  to  study  the  0(  P)  +  HF 
and  0(  P)  +  HC1  collisions.  These  should  have  most  of  the  same  properties 
already  observed  for  H2O  and  CC^,  but  may  add  some  new  features  since  HOF, 
HFO,  etc.,  probably  have  a  degree  of  stability.  The  PES  for  these 
surfaces  are  easier  to  obtain  precisely  because  they  have  only  three 
degrees  of  freedom,  compared  to  six  for  the  H2O  and  CO2  collisions. 

Again,  the  experiments  are  highly  difficult  for  these  systems,  and 
reliable  calculations  should  be  a  cost  effective  approach  to  obtaining 
the  necessary  kinetic  data. 

A  fourth  direction  for  additional  work  would  be  the  H2O  with 
HjO  collision.  This  seems  to  be  an  important  contributor  to  plume 
signatures.  Since  it  involves  12  degrees  of  freedom,  and  one  can  compute 
only  a  limited  number  of  points  on  the  surface,  the  behavior  of  the 
analytic  fit.  to  this  surface  would  have  to  account  for  much  more  of 
the  shape  of  the  surface  than  in  the  0  +  H2O  or  COp  collisions.  However, 
we  have  an  exceptional  fit  for  H2O  itself  which  could  be  incorporated 
into  the  six-atom  surface.  In  addition,  there  is  only  a  single  PES  for 
H2O  +  H2O  and  no  open-shells  are  involved  to  complicate  the  calculations. 
The  H2O  +  H2O  surface  would  have  a  number  of  local  minima  that  the  fit 
would  need  to  correctly  describe. 


One  other  problem  worth  mentioning  is  the  collision  of 
with  H2O  and  CO^ -  Again,  there  is  only  a  single  PES  for  their  closed- 
shell  molecules,  and  there  is  quite  a  bit  of  experimental  information 
on  this  system.  It  would  be  Interesting  to  study  these  collisions  as  a 
function  of  vibrational/rotational  state  for  N^,. 

In  the  area  of  plume  identification  and  detection,  many  other 
possibilities  for  future  theoretical  calculations  exist,  since  much  of 
the  resolution  of  the  signature  problem  rests  upon  a  detailed  knowledge 
of  the  atom-atom,  atom-molecule,  and  molecule-molecule  interactions. 
First  principle  theoretical  chemistry  can  contribute  significantly  to 
all  of  these  problems . 
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